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Executive  Summary 


The  integral  equation  constrained  optimization  approach  to  finding  three-dimensional  minimum-drag  shapes 
for  bodies  translating  in  viscous  incompressible  fluid  has  been  developed.  The  approach  relies  on  the  theory  of 
generalized  analytic  functions  for  obtaining  efficient  integral  equations  for  3D  boundary-value  problems  with 
linear-  partial  differential  equations  (PDEs),  e.g.,  the  Stokes  and  Oseen  equations  for  viscous  incompressible 
fluid,  linearized  equations  of  magnetohydrodynamics  (MHD)  etc.  It  assumes  the  following  steps:  (i)  identify¬ 
ing  a  class  of  generalized  analytic  functions  related  to  the  PDEs;  (ii)  representing  corresponding  fields,  e.g., 
the  fluid  velocity  field,  the  electromagnetic  field  etc.,  in  terms  of  generalized  analytic  functions  from  the  identi¬ 
fied  class;  and  (iii)  reformulating  boundary-value  problems  for  the  involved  fields  as  boundary- value  problems 
for  generalized  analytic  functions,  which  could  then  be  reduced  to  integral  equations  based  on  the  generalized 
Cauchy  integral  formula.  Solutions  of  the  obtained  integral  equations  can  be  represented  by  finite  function 
series  with  series  coefficients  determined  by  quadratic  error  minimization.  The  framework  of  generalized  an¬ 
alytic  functions  has  been  developed  in  application  to  3D  Stokes  and  Oseen  flows,  two-phase  Stokes  flows  and 
3D  magnetohydrodynamic  flows  governed  by  linearized  MHD  equations. 

In  shape  optimization  problems,  the  suggested  approach  replaces  3D  boundary-value  problems  with  govern¬ 
ing  PDEs  by  corresponding  boundary  integral  equations.  Minimum-drag  shapes,  represented  in  finite  function 
series  form,  are  then  found  by  the  adjoint  equation-based  method  with  a  gradient-based  algorithm,  in  which 
the  gradient  for  shape  series  coefficients  is  determined  analytically.  Compared  to  PDE  constrained  optimiza¬ 
tion  coupled  with  the  finite  element  method  (FEM),  the  approach  reduces  dimensionality  of  the  flow  problems, 
solves  the  issue  with  region  truncation  in  exterior  problems,  finds  minimum-drag  shapes  in  semi-analytical 
form,  and  has  fast  convergence.  Its  efficiency  has  been  demonstrated  in  solving  three  drag  minimization  prob¬ 
lems  under  the  Oseen  approximation  of  the  Navier-Stokes  equations  for  different  Reynolds  numbers:  (i)  for  a 
body  of  constant  volume,  (ii)  for  a  torpedo  with  only  fore  and  aft  noses  being  optimized,  and  (iii)  for  a  body 
of  constant  volume  following  another  body  of  fixed  shape  (e.g.  torpedo  chasing  a  target).  The  minimum-drag 
shapes  in  problem  (i)  are  in  good  agreement  with  the  existing  optimality  conditions  and  conform  to  those  ob¬ 
tained  by  PDE  constrained  optimization.  Problem  (ii)  has  shown  that  the  minimum-drag  shape  for  the  torpedo 
is  fore-aft-symmetric,  whereas  problem  (iii)  has  revealed  that  the  minimum-drag  shape  for  the  trailing  body  is 
only  slightly  more  prolate  than  the  one  in  problem  (i)  for  same  Reynolds  number. 

For  MHD  flows  in  the  presence  of  aligned  magnetic  fields,  the  necessary  optimality  condition  for  the 
minimum-drag  shapes  subject  to  a  volume  constraint  has  been  obtained  analytically.  It  has  been  shown  that  re¬ 
gardless  of  magnitudes  of  Hartmann  number,  Reynolds  number,  and  magnetic  Reynolds  number,  the  minimum- 
drag  shapes  are  fore-and-aft  symmetric  and  have  conic  endpoints  with  the  angle  of  2ji/3.  For  fixed  Reynolds 
and  magnetic  Reynolds  numbers,  the  drag  reduction  as  a  function  of  the  Cowling  number  is  smallest  at  1.  In 
considered  MHD  problems,  the  drag  coefficients  for  the  minimum-drag  shapes  and  minimum-drag  spheroids 
are  sufficiently  close. 

The  project  involved  Postdoctoral  Associate  Anton  Molyboha  from  the  Department  of  Mathematical  Sci¬ 
ences,  Stevens  Institute  of  Technology. 
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Shape  Optimization  in  Magnetohydrodynamics 
1  Introduction 

Industrial  applications  of  magnetohydrodynamics  (MHD)  include  but  are  not  limited  to  MHD  flow  control 
schemes  for  hypersonic  vehicles  |fl9l,  torpedo  drag  reduction  with  MHD  boundary-layer  control  |[2l,  and  MHD 
propulsion  systems  lf25l. 

The  steady  flow  of  an  electrically  conducting  viscous  incompressible  fluid  in  the  presence  of  magnetic  field 
and  with  neglected  thermal  effects  can  be  characterized  by  three  independent  parameters:  Hartmann  number 
M,  Reynolds  number  R,  and  magnetic  Reynolds  number  Rm.  The  Hartmann  number  is  the  ratio  of  the  Lorentz 
force  to  the  viscous  force  in  the  Navier-Stokes  equations  (when  M  =  0,  the  velocity  field  is  uncoupled  from  the 
electromagnetic  field),  whereas  the  magnetic  Reynolds  number  is  interpreted  as  the  ratio  of  magnetic  advection 
to  magnetic  diffusion  in  the  combination  of  Ohm’s  and  Ampere’s  laws  (when  Rm  =  0,  the  magnetic  field  is 
uncoupled  from  the  velocity  field).  The  ratio  of  the  magnetic  forces  to  the  inertial  forces  is  characterized  by  the 
Cowling  number  sQwhich  can  be  expressed  as  M2/ (Rm  R)  when  R  /  0  and  Rm  ^  0. 

An  MHD  problem  that  has  attracted  much  of  the  attention  is  arguably  the  one  of  an  electrically  conducting 
flow  past  a  nonmagnetic  sphere  in  the  presence  of  a  uniform  magnetic  field  being  aligned  with  the  undisturbed 
flow;  see  e.g.  H  Q21I271  ED  [01.  For  this  problem.  Figure  [T]  shows  the  drag  for  the  sphere  normalized  to  the 
sphere  Stokes  drag  as  a  function  of  S  for  R  =  Rm  =  1  (curve  a);  R  =  1,  Rm  =  3  (curve  b);  R  =  Rm  =  2  (curve 
c);  R  =  3,  Rm  =  1  (curve  d);  and  R  =  Rm  =  3  (curve  e).  At  S  =  1,  the  drag  attains  minimum  and  is  nonsmooth; 
see  MM.  Namely  the  fact  that  the  sphere  drag  at  S  =  1  is  nonsmooth  is  remarkable.  It  poses  the  research 
questions:  If  the  sphere  is  replaced  by  an  arbitrary  nonmagnetic  body  of  revolution,  what  is  body’s  shape  that 
has  the  smallest  drag  subject  to  a  volume  constraint  and  how  does  it  depend  on  S?  Are  there  any  qualitative 
differences  in  the  minimum-drag  shapes  for  S  <  1,  S  =  1,  and  S  >  1?  How  does  drag  reduction  depend  on  S? 
Answering  these  questions  is  the  subject  of  this  study. 


Figure  1 :  Drag  for  the  unit  sphere  normalized  to  the  sphere  Stokes  drag  as  a  function  of  S  for  R  =  Rm  =  1 
(curve  a);  R  =  1,  Rm  =  3  (curve  b);  R  =  Rm  =  2  (curve  c);  R  =  3,  Rm  =  1  (curve  d);  and  R  =  Rm  =  3  (curve  e). 

The  challenge  in  addressing  the  posed  questions  is  that  the  traditional  PDE  constrained  optimization  ap¬ 
proach  coupled  with  the  finite  element  method  (FEM)  is  slowly  converging  and  inaccurate  in  general.  The 
second  deficiency  is  attributed  to  the  fact  that  being  applied  to  external  problems,  the  approach  truncates  and 
discretizes  an  external  domain  and  finds  minimum-drag  shapes  point  wisely. 

Our  approach  reduces  the  MHD  problem  to  boundary  integral  equations,  derives  the  optimality  condition 
for  the  minimum-drag  shapes  analytically,  and  obtains  minimum-drag  shapes  in  a  functional  series  form.  To 

'This  number  is  also  denoted  by  Co  and  is  also  called  “pressure”  number  in  127 1 1 13J. 
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this  end,  the  MHD  equations  are  linearized,  and  the  velocity,  pressure  and  magnetic  fields  in  and  out  the 
body  are  represented  by  four  generalized  analytic  functions.  Under  the  assumption  that  the  fluid  and  body 
are  both  nonmagnetic  and  have  same  magnetic  permeability,  the  axially  symmetric  MHD  problem  is  reduced 
to  integral  equations  for  the  boundary  values  of  two  generalized  analytic  functions  based  on  the  generalized 
Cauchy  integral  formula.  The  Hartmann  number  is  assumed  to  be  nonzero  since  when  M  =  0,  the  velocity  field 
is  uncoupled  from  the  magnetic  field.  With  M  /  0,  three  cases  are  studied  separately:  (a)  Rm  /  0,  RmR  /  M2 
(S  /  1);  (b)  Rm  =  0;  and  (c)  RmR  =  M2  (S  =  1).  The  reason  for  this  is  that  in  each  case,  solving  the  MHD 
problem  involves  different  number  of  generalized  analytic  functions  from  two  classes:  r- analytic  functions  and 
//-analytic  functions.  This  is  the  mathematical  explanation  of  why  the  case  S  =  1  is  special. 


1.1  r-Analytic  and  //-Analytic  Functions 


A  function  F(x,y )  =  U  (x,y)  +i  V(x,y)  with  i  =  y/—l  is  called  generalized  analytic  or  pseudo-analytic  if  its  real 
and  imaginary  parts,  U  and  V,  respectively,  satisfy  the  so-called  Carleman  or  Bers-Vekua  system  |4]  [26] 


dU 

dx 


-=r — \-aU  +  bV  =  0. 
dy 


dU  dV  , 

- b  — VcU  -\~dV  —  0, 

dy  ax 


(1) 


where  a  =  a(x.y),  b  =  b(x,y),  c  =  c(x,y),  and  d  =  d(x,y)  are  real-valued  functions.  For  example,  for  a  = 
b  =  c  =  d  =  0,  the  system  0  defines  ordinary  analytic  functions.  Generalized  analytic  functions  arise  in 
various  areas  of  applied  mathematics  including  hydrodynamics,  gas  dynamics,  theory  of  elasticity,  heat  transfer, 
electromagnetism,  quantum  mechanics,  etc.;  see  J61I26,  2J.I  TJ .  Their  theory  has  been  extensively  developed 
since  the  mid-20th  century  and  extends  the  majority  of  results  for  ordinary  analytic  functions,  e.g.  the  formal 
powers  |5]l!l  and  the  Cauchy  integral  formula  li4ll251|2Hf8l. 

Several  important  classes  of  generalized  analytic  functions  arise  from  the  relationship 


curl  $+[ax$]  =  -  gradT1,  div  $  =  0 


(2) 


for  a  vector  field  $  and  scalar  field  VF,  where  a  is  a  known  vector  function.  This  relationship  is  frequently 
encountered  in  problems  of  applied  mathematics.  For  example,  for  T1  =  0  and  a  =  0,  <[2])  simplifies  to  an 
irrotational  solenoidal  field  $  found  in  electrostatics  and  ideal  fluid,  whereas  for  a  =  0,  it  defines  so-called 
related  potentials  lF  and  <F,  e.g.  the  pressure  and  vorticity  in  the  Stokes  equations,  electric  potential  and 
magnetic  field  in  conductive  materials,  etc.;  see  EE  For  a  y/0  and  T*  7/ 0,Q  arises  in  the  Maxwell  equations 
for  quasi-stationary  electromagnetic  fields,  in  the  Oseen  equations  and  linearized  MHD  equations;  see  ll30l 
Example  4]. 

Let  (r,  tp, z)  be  a  cylindrical  coordinate  system  with  the  basis  (er,e<p,k).  In  the  axially  symmetric  case  with 
the  z  axis  being  the  axis  of  revolution,  let  $  =  Ty(r.z)  e,.  +  d>-{r.z)  k.  VP  =  0,  and  a  =  — 2Ak,  where  A,  is  a 
constant.  Then  for  U  =  e  and  V  =  e  kz't>l,  0  reduces  to 


dU 

dr 


+  A  [/ 


which  is  a  special  case  of  0  and  implies  that 

( Aq  —  A2)  U  =  0  and 


(a1-a2)f  =  0, 


(3) 

(4) 


where  A*  denotes  the  so-called  k-harmonic  operator:  A*  =  d1  /dr2  +  r'1  d/dr  +  d2/d~2  —  k2/r2.  The  functions 
U  and  V  satisfying  0  form  an  //-analytic  function  G  =  U  +  iV  of  a  complex  variable  C,  =  r  +  \z,.  For  A  =  0, 
0  defines  a  so-called  r-analytic  function  and  is  arguably  the  most  studied  system  among  various  classes  of 
generalized  analytic  functions.  It  arises  in  the  axially  symmetric  theory  of  elasticity  GT1  [Q  and  in  the  axially 
symmetric  Stokes  and  Oseen  flows  ll35l[34i  i3 1 , 28  29l  301 .  Since  an  r-analytic  function  satisfies  0  for  A  =  0,  it 
is  also  referred  to  as  0- harmonically  analytic  function  |[28ll.  Both  classes  of  r-analytic  and  //-analytic  functions 
are  instrumental  in  constructing  solutions  to  axially  symmetric  MHD  problems. 
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1.2  Generalized  Cauchy  Integral  Formula  and  Series  Representation 

Let  G+  be  a  generalized  analytic  function  in  a  bounded  open  region  D 1  in  the  right-half  /"-plane  (D '  may 
contain  parts  of  the  z  axis).  The  boundary  of  D  is  assumed  to  be  a  piece-wise  smooth  positively  oriented 
curve  l,  which  is  either  closed  or  open  with  the  endpoints  lying  on  the  z  axisj^J  Let  D  be  the  complement  of 
D+u£  in  the  right-half  rz- plane  CD  is  unbounded),  and  let  G  be  a  generalized  analytic  function  in  D  that 
vanishes  at  infinity.  For  convenience,  an  arbitrary  function  f(r,z)  will  be  denoted  by  f(Q  without  assuming 
its  analyticity.  Let  G±(Q  satisfy  the  Holder  condi tior0 on  /,  and  let  If  be  the  reflection  of  i  over  the  z  axis. 
Then  under  the  assumption  of  the  symmetry  condition  G±(— Q  =  G±(Q,  the  Cauchy  integral  formula  for  G± 
is  given  by  ITI  j2T.  261: 

G±{Q  =  ±^-<b  G±(z)W(^)dz,  C,  G  D±,  (5) 

Jh \jg 

where  £  =  r+iz,  x  =  r\  +  izi  and  TF(^,x)  =  W(r,z\ri,zi)  is  a  generalized  Cauchy  kernel.  If  the  the  boundary 
l  in  smooth,  then  on  £,  G±(Q  satisfies  the  generalized  Sokhotski-Plemelj  formula: 


G±(Q  =  V(Q  G±(T)TF(^x)dx,  C  G  *■ 

2  2n\Jt\jii 


(6) 


If  £  is  piece-wise  smooth  (has  salient  points),  in  particular  if  the  endpoints  of  l  lie  on  the  z  axis  and  the  angle 
between  £  and  the  z  axis  at  one  of  the  endpoints  is  not  tc/2  (conic  endpoint),  then  at  a  salient  point  of  l,  the 
coefficient  1/2  in  the  generalized  Sokhotski-Plemelj  forlmula  (|6])  is  replaced  by  a  finite  function  of  the  angle 
between  tangents  at  the  salient  point;  see  [QJ  (31.13a),  (31.13b)]  and  ll3~31  (3.7)]. 

Let  Wr  and  Wh  denote  the  generalized  Cauchy  kernels  for  r-analytic  and  //-analytic  functions,  respectively. 
For  r-analytic  functions,  ‘Wr  is  given  in  HI,  ll28l  Theorem  2]  and  ll33l  Theorem  3.2].  Theorem  1  and  Corollary 
4  in  lf30l  present  two  forms  for  ‘Wu  with  X  >  0.  The  next  result  extends  the  Cauchy  integral  formula  for 
H- analytic  functions  for  arbitrary  X  f  0. 


Theorem  1  Let  G±  be  an  H -analytic  function  in  D+  that  satisfies  Q  with  X  f  0  and  satisfies  the  Holder 
condition  on  t.  Suppose  G±(— Q  =  G±(Q.  Then  the  generalized  Cauchy  integral  formula  ©  takes  the  form 


GH  q 

where 


1 

2jti 


G(x)TFH(£,x)dx 


f Kl^1x,X)Re[G±(x)dz}+iK2(^x,X)lm[G±(x)dx}  ,  ^eD±, 
zm  j  t 

(V) 


Ki(W 0  =  -n 


r *  e-|k|p(C,T,0 

L  p(C,v) 


( (C  -  n  cosf  +  izi)  (1  +  |X|  p(C,  V)) 

^  p(?,V)2 3 


dt, 


K-2&,  T,X)  =  -n 


e  .  (zi-z)cosf 

- (/-cost  —  ri)  (1  +  |A,|p(£,x,t))  +t- - 


p(C.x,0: 


e“|k|p(£>x>f)  _  e^(zi-*) 

+  i - —Til - r -  I  COSt  + 


A(C,x,t) 


P(C,x,02 

(r  — n  cosf)(n  —  rcosf)(2A(£,x,t) -z\  +z) 
p(^,x,t)2A(^x,t) 


(/•—  /'I  COSt)(/-j  —  /-cost) 


-Fie  (2A,cost  +  |A,|- 

V  P(C,x,f)A(C,x,t) 


dt 


J  P(^V)' 


p(C,X,t)  =  \Jr2  +  rf  -  2/ti  cost  +  (z  -  zi  )2, 
A(?,x,t)  =zi  -z  +  p(^,x,t)sign^. 


2Open  segments  of  the  z  axis  that  ©+  may  contain  are  not  parts  of  £. 

3This  condition  means  that  for  some  parametrization  ^(/)  of  l,  the  boundary  value  f{C,(t))  satisfies  |/(C(h)) — /(C(h))l  <  c |fi  — t\ 
for  all  t\  and  t2,  some  p  e  (0, 1],  and  nonnegative  constant  c. 
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Proof.  The  proof  is  similar  to  that  of  Theorem  1  in  iPOl.  □ 

For  the  region  exterior  to  a  sphere.  Example  5  and  Proposition  2  in  11301  provide  series  representations  for 
r-analytic  functions  and  //-analytic  functions  with  X  >  0,  respectively.  The  next  proposition  generalizes  the 
representation  for  //-analytic  functions  for  any  X  f  0. 

Proposition  1  (//-analytic  function  in  the  region  exterior  to  a  sphere)  For  the  region  exterior  to  a  sphere 
centered  at  the  origin,  an  H -analytic  function  satisfying  ([3])  for  Xf^O  and  vanishing  at  infinity  can  be  rep¬ 
resented  in  the  spherical  coordinates  (A\  th  cp)  by 

H{R, ■&)  =R~^Y2=\An  (l„(cosi3)K,i+i(|X|/?)  -sign/iL_„(cosO)Kn_i(|)i|/?))  ,  (8) 

where  L„  (cos  tt)  =  n  P„  (cos  tt)  —  i  P„  (cos  0),  P„  (cos  if)  is  the  associated  Legendre  polynomial  of  the  first  kind 
of  order  n  and  rank  k  ( for  k  =  0,  the  superscript  is  omitted ),  KJ1+ 1  (•)  is  a  modified  spherical  Bessel  function  of 
the  third  kind  (see  Sec.  7.2.6]),  and  real-valued  coefficients  An  are  such  that  the  corresponding  series  © 
converge  for  all  0  G  [0, 7C]  and  R  greater  than  or  equal  to  the  radius  of  the  sphere. 

Proof.  The  proof  is  analogous  to  the  proof  of  Proposition  2  in  Ii30ll .  □ 


2  Magnetohydrodynamics 

The  steady  flow  of  an  electrically  conducting  viscous  incompressible  fluid  is  governed  by  the  MHD  equations 

'  p(U-grad)U  =  —  gradp  +  pvAU+ju[J  xH], 

<  curlH  =  471 J,  J  =  c(E+q[UxH]),  (9) 

curlE  =  0,  divU  =  0,  divH  =  0, 

where  U  is  the  fluid  velocity,  p  is  the  pressure  in  the  fluid,  E  and  H  are  the  electric  and  magnetic  fields, 
respectively,  J  is  the  current  density,  v  is  the  kinematic  viscosity,  p  is  the  fluid  density,  p  is  the  magnetic 
permeability,  and  a  is  the  conductivity;  see  l22ll23:T61. 

Suppose  a  nonmagnetic  solid  body  of  revolution  is  immersed  in  the  uniform  flow  in  the  presence  of  mag¬ 
netic  field.  At  infinity,  the  flow  and  magnetic  field  are  assumed  to  be  constant  and  aligned  with  the  z-axis,  i.e. 
UU  =  Took  and  Hj^  =  //<*, k.  where  VT  and  //<*,  arc  constants.  Also,  the  body’s  axis  of  revolution  is  parallel  to 
the  z-axis.  Under  different  assumptions,  this  MHD  problem  was  considered  in  ll9l  IT3l IT4l  fl~5l  fTTl l24l ITOl .  Let  u 
be  the  disturbance  of  the  velocity  field:  U  =  VU(k  +  u),  and  let  h+  and  h  be  the  disturbances  of  the  magnetic 
field  in  and  out  the  body,  respectively:  H  =  H«,  (k  +  h+)  in  the  body  and  H  =  H X  (k  +  h~ )  out  the  body.  On  the 
boundary  S  of  the  body  and  at  infinity,  u  and  p  satisfy  the  conditions 

uls  =  — k,  u|oo  =  0,  PL  =  0.  (10) 

Inside  the  body,  the  electric  field  is  zero,  whereas  the  magnetic  field  satisfies  curlH  =  0  and  divH  =  0,  or 
in  terms  of  the  disturbance  h+, 

curlh+  =  0,  divh+  =  0.  (11) 

As  in  it  is  assumed  that  the  fluid  and  body  are  nonmagnetic  and  have  same  magnetic  permeability.  In 

this  case,  the  magnetic  field  is  continuous  across  the  boundary  of  the  bod\J^]and  also  h  vanishes  at  infinity: 

h+L  =  h1s’  h"L  =  0.  (12) 

4If  the  fluid  and  body  have  different  magnetic  permeability  then  across  the  body’s  boundary,  the  tangential  component  of  the 
magnetic  field  H  and  the  normal  component  of  the  magnetic  induction  B  should  be  continuous. 
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The  problem  (|9l)-(fT2|)  is  axially  symmetric.  In  this  case,  the  velocity,  pressure  and  electromagnetic  field  are 
independent  of  the  angular  coordinate  cp,  i.e. 

u  =  ur(r,z)er  +  uz(r,z)K  p=p(r,z),  h±  =  hf(r,z)er  +  hf(r,z)k,  E  =  £,(p(r,z)e(p. 

The  equation  curlE  =  0  implies  that  rE <p  =  c,  where  c  is  a  constant.  Since  [u  x  H]  and  curlH  vanish  at  infinity, 
we  have  c  =  0,  and  consequently,  E  =  0. 

Rescaling  the  linear  dimensions,  velocity,  pressure,  and  magnetic  field  by  a,  14,,  V4pv/u,  and  Hoo,  respec¬ 
tively,  where  a  is  the  half  of  the  diameter  of  the  body,  and  assuming  u  and  h  to  be  small,  we  can  rewrite 
equations  ((9j)  without  E  and  J  in  the  dimensionless  linearized  form 

’  R  (k  -grad)u  =  -grad£?+  Au  +  M2[[(u  —  h~)  x  k]  x  k], 

<  curllG  =  Rm  [(u  —  h~)  x  k],  (13) 

divu  =  0.  divh~  =  0. 

where  R  =  V4  fl/v  is  the  Reynolds  number,  Rm  =  47tV4  ape  is  the  magnetic  Reynolds  number,  and  M  = 
jj  //<*,  a v/a/(pv)  is  the  Hartmann  number. 

For  M  =  0,  the  first  equation  in  (13])  reduces  to  the  Oseen  equations,  and  the  problem  for  the  velocity  and 
pressure  becomes  uncoupled  from  the  magnetic  field.  For  solving  the  Oseen  equations,  see  e.g.  Ii30ll.  For 
M  /  0,  there  are  three  cases  to  analyze: 

(a)  Rm  /  0  and  RmR  /  M2:  A  solution  to  (fT3|)  is  represented  by  two  //-analytic  functions  and  one  r-analytic 
function,  whereas  a  solution  to  (fTT|)  is  given  by  a  single  r-analytic  function. 

(b)  Rm  =  0:  The  magnetic  field  is  constant  everywhere,  i.e.  h+  =  0,  and  a  solution  to  (f]~3j)  is  represented  by 
two  //-analytic  functions.  (The  case  of  R  =  0  yields  no  simplification.) 

(c)  RmR  =  M2  (S  =  1):  A  solution  to  (|T3j)  is  represented  by  one  //-analytic  function  and  two  r-analytic 
functions,  and  a  solution  to  (fTTj)  is  given  by  a  single  r-analytic  function. 

Let  £  be  the  positively  oriented  cross  section  of  the  body  in  the  right-half  rz-plane  (r  >  0),  and  let  £'  be  the 
reflection  of  £  over  the  z  axis.  The  curve  £  divides  the  right-half  rz-plane  into  interior  bounded  region  ‘D 1  and 
exterior  unbounded  region  ‘D  . 

Theorem  2  (solution  representation,  case  (a))  In  the  axially  symmetric  case  with  M  /  0  and  RmR  /  M2,  a 
solution  to  <UU)  and  CD  is  given  by 

UZ  +  i ur  =  ((1-2X2*)  ex* 1  Gi  -  ( 1  -  2Xl  *)  G2)  -  R  R*_  Gf ,  (14) 

K-+ih'  =  (^Gl -  R <15> 

h+  +  ih+  =  G+,  (16) 

where  G\  and  G2  are  H -analytic  functions  in  ‘D  ~  satisfying  ([3])  with  k  =  A,|  2  =  R+Rm±^R  Rm^  +1D(L1  respec¬ 
tively,  and  vanishing  at  infinity;  Gt  and  Gf  are  r-analytic  functions  in  ‘D  an  D~ ,  respectively,  with  G 
vanishing  at  infinity;  and  x  =  Rm/  (RmR  —  M2). 

In  this  case,  the  pressure  and  scalar  vortex  function  or  =  e(p  •  curlu  (i.e.  curlu  =  ore<p)  are  determined  by 
P=  (^e^Re^]  -Xie^R e[G2])  +  Re[G3"],  (17) 

co  =  ((R-2^)eXlZIm[G1]-(R-2Xi)e^2Im[G2])  .  (18) 
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Proof.  With  the  identity  (k  •  grad)  u  =  —  curl  [k  x  u],  the  first  equation  in  (fT3]l  can  be  rewritten  in  the  form 

curl  (to  —  R[k  x  u])  —  M2  [k  x  [k  x  (u  —  )]]  +  grad£?  =  0,  (19) 

where  to  =  curlu  is  the  vorticity.  The  continuity  equation  divu  =  0  implies  that  u  and  to  satisfy  the  identity 

curlfk  x  u]  +  grad(k  •  u)  —  [k  x  to]  =  0,  (20) 

whereas  div  h  =  0  implies  a  similar  identity  for  h 

[k  x  curlh~]  =  curl[k  x  h_]  +  grad(k  •  h“), 

with  which  the  second  equation  in  (p~3|)  can  be  represented  in  the  form: 

curlfk  x  h“]  +  grad(k-  h  )  +  Rm  [k  x  [k  x  (u  —  h-)]]  =  0.  (21) 

Forming  a  linear  combination  of  (|T9|),  (|20|),  and  (|2T])  with  constant  weights  pi,  P2,  and  P3,  respectively,  we  have 

curl(Pito  +  [k  x  ((p2  -  Pi  R)  u  +  p3 h“)])  -  [k  x  (p2to  +  (Pi  M2  -  p3Rm)[k  x  (u  -  h")])] 

+  grad(pip+p2(k-u)  +  p3(k-h"))  =  0. 

^  Let  p3  =  -(p2  -  Pi  R)  and  (p2  -  pi  R)/px  =  (pi  M2  -  p3Rm)/p2,  then  (P2/P1)2  -  (R  +  Rm)p2/Pt  +  RmR  - 
M2  =  0,  whence  P2/P1  =  2 A*,  k  =  1,2,  and  ([22])  becomes  a  particular  case  of  ([2]): 

curlAit  —  2A,;fc[k  x  Ax] +grad'Pjt  =  0,  divA/;  =  0,  k=  1,2,  (23) 

where 

A/t  =  to  +  (2A*-R)[kx(u-h-)]J  'P*  =  p+2A*(k-u)  +  (R-2A*)(k-h-).  (24) 

In  the  axially  symmetric  case,  ([23])  reduces  to  the  system  ([3])  determining  //-analytic  functions.  Let  A/c  =  A/{  e<p. 
k  =  1,2,  and  to  =  coe<p.  In  this  case, 

'¥k  +  iAk  =  2eXkZGk,  k  =  1,2,  (25) 


where  Gk  are  //-analytic  functions  defined  in  the  theorem. 

On  the  other  hand,  if  p2  =  0  and  Pi  M2  —  p3  Rm  =  0  then  ([22])  becomes  the  relationship  for  related  potentials'. 

curl  A3  +  grad  'Pj  =  0,  div  A3  =  0,  (26) 


where 

A3  =  Rmto  +  [k  x  (M2h~  -  RmRu)],  T/3  =  Rmp+  M2(k-h_). 

In  the  axially  symmetric  case,  A3  =  A3  e(f)  and  ([26])  determines  an  /  -analytic  function  C3  in  ‘D  : 


(27) 


^3  +  1^3  =  G, 


Equations  (23 )— ( 28 )  can  be  reformulated  in  the  complex  form 


(28) 


2e^'  Gk  —  p  +  i(0  +  2 -f* in,-)  —  i R ur  T  (R  —  2A,yt) (/r~  T ihr  ) ,  k  —  1,2, 

G3  =  Rm(^?+ iro)  —  i R Rm ur 4-  M"(/i_  +i/?r), 

from  which  the  representations  (fT4[),  (fi3]l,  (17]),  and  (fl8[)  follow.  The  condition  M  /  0  guarantees  that  \\  /  L2. 

Finally,  in  the  axially  symmetric  case,  h+,  being  an  irrotational  solenoidal  field  (see  CD).  can  be  repre¬ 
sented  by  (fl6[).  □ 


Corollary  1  (solution  representation:  case  (b))  In  the  axially  symmetric  case  with  M  /  0  and  Rm  =  0,  the 
magnetic  disturbance  is  zero,  i.e.  h+  =  0,  and  the  velocity,  pressure,  and  vorticity  that  solve  ©  can  be 
represented  by 


uz  +  \ur  =  (e^zGi  -e^zG2)  , 

»,+i“>  =  el,;(G'-2(XGl5B|)+el“ 


R 

2(Xi  —  A2) 


(29) 


where  G\  and  G2  are  H -analytic  functions  satisfying  ([3])  with  X,  =  Ai>2  =  ^  /R±  \/R2  +  4M2j,  respectively, 
and  vanishing  at  infinity. 


Detail.  For  Rm  =  0,  (jT3]>  reduces  to  h7  +i hf  =  M~2  Gf  which  with  (|T6j)  and  (JT2])  implies  M~2 Gf  =  G 3"  on  l. 
Consequently,  since  Gf  vanishes  at  infinity,  by  the  Sokhotski-Plemelj  formula  Q,  Gf  =  0  in  D±,  respectively, 
and  (f2^>  follows  from  (fT4(),  (fTT]),  and  ([T8l).  □ 


Corollary  2  For  M  f  0  and  Rm  =  R  =  0  the  representation  ( 29 1  simplifies  to 

uz  +  iur  =  ^  ^eMz/2Gi  _ e_Mz/2G2j  ,  £?+ico  =  eMz/2Gi  +e_MZ//2G2, 

where  G\  and  G?  are  H -analytic  functions  satisfying  ([3])  with  X  =  M/2  and  X  =  —  M/2,  respectively.  It  is 
similar  to  the  solution  form  suggested  by  Chester  /[9]/  in  the  case  of  R  =  Rm  =  0. 


Corollary  3  For  M  =  0  and  R  f  0,  the  representation  (29)  reduces  to  the  representation  (22)— (23)  in  l\30\l  for 
the  velocity,  pressure  and  vorticity  for  the  axially  symmetric  Oseen  flow  of  nonconducting  fluid. 

Theorem  3  (solution  representation,  case  (c))  In  the  axially  symmetric  case  with  M  fl  0  and  Rm  R  =  M2,  a 
solution  to  ©  and  ©  is  given  by 


lly  +  i  ur  = 


1 


2R 


R  +  Rm  V  R  +  R 


Rn 


e^zGi+  (  R 


R 


Z  2/  R+Rn 


h-.+'K  =  R+K 


R+Rn 


eAzGi  +  I  z--r+ 


1 


R+Rn 


G2  +  RmG3  1  , 

G2  +  G3), 


lC  +  ihf  =  Gf, 


(30) 


(31) 


(32) 


where  X  =  (R  +  Rm)/2,  G\  is  an  H -analytic  function  in  ‘D  that  satisfies  ([3])  with  X  and  vanishes  at  infinity;  G2 
and  Gf  are  r -analytic  functions  in  ‘D  and  vanishing  at  infinity;  and  G/  is  an  r-analytic  function  in  F)+. 

In  this  case,  the  pressure  and  vortex  function  are  determined  by 


P  = 


RmR 


R  +  Rm  V  R  +  R 


eXzRe[Gi]-Re 


z-r+RTi^ 


g2 


—  Re [G3  ]  )  +  Re[G2] , 


to  =  (2ReXzIm[Gi]  +  Rm  Im[G2])  . 


(33) 

(34) 


Proof.  The  proof  is  partially  based  on  the  proof  of  Theorem^  When  RmR  =  M2,  the  relationships 
hold  for  An  =  (R  +  Rm)/2  and  Xi  =  0.  For  An  =  (R  +  Rm)/2,  ((251)  simplifies  to 


-(25) 


2e(R+Rm )-/- (j j  —  +- i co  +  ( R  +  Rm) (mz  +  in,.)  —  i Rur  —  Rm (hz  +  ihr  ), 


(35) 
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where  G i  is  the  //-analytic  function  defined  in  this  theorem.  The  case  of  L2  =  0  means  that  (|23[)  reduces  to  the 
relationship  (|26|)  for  the  related  potentials.  Consequently,  a  different  approach  is  required.  Since  RmR  =  M2, 
multiplying  (|2T|)  by  R  and  adding  to  (p~9|).  we  obtain 

curl(cj  -  R[k  x  (u  — h~)])  +  grad(£?+  R(k-  h-))  =  0, 
which  with  the  second  equation  in  ([13])  and  conditions  div  u  =  0  and  div  h  =  0  can  be  rewritten  as 


A(u  +  (R/Rm)h  )  =  gradO+R(k-h  )),  div(u+ (R/Rm)h  )  =  0. 


(36) 


The  system  (f36[)  is  similar  to  the  Stokes  equations  for  the  viscous  incompressible  fluid.  In  the  axially 
symmetric  case,  its  solution  is  given  by  Proposition  7  in 


R 


uz  +  iur+—(hz  +i  hr)  =  [z-  -r  )  G2  +  G3  , 


ico  —  i R iiy f-  R(/r  T xhr )  —  G2 , 


(37) 


where  G2  and  G3  are  the  r-analytic  functions  defined  in  this  theorem.  Observe  that  the  conditions  M  /  0  and 
RmR  =  M2  guarantee  that  Rm  /  0.  The  representations  (j30j),  (pT|),  (|33|),  and  ([34])  follow  from  ([35])  and  (|37|), 
whereas  ([32])  remains  same  as  in  Theorem [2]  □ 

The  advantage  of  the  representations  ([T4|)-([T8]),  (29),  and  (30 (-([34])  compared  to  the  existing  forms,  e.g. 
in  @,  is  that  these  representations,  being  linear  combinations  of  the  generalized  analytic  functions,  involve 
no  derivates  of  unknown  functions  and  represent  simultaneously  all  three:  the  velocity,  vorticity  and  pressure. 
This  fact  considerably  simplifies  deriving  boundary  integral  equations  based  on  the  generalized  Cauchy  integral 
formula. 


3  Boundary-value  Problems  and  Integral  Equations 

3.1  Case  (a):  M^O,  Rm^0,  and  RmR  ^  M2 

For  M  /  0  and  RmR  /  M2,  equations  (T3|)  with  the  boundary  conditions  ([TO])  and  (12)  reduce  to  the  boundary- 
value  problem  for  the  four  generalized  analytic  functions  G\,  G2,  and  Gf  defined  in  Theorem|2j 


(1  -2X2x)e^K  Gi  —  (1  —  2Xix)e^2'  G2  —  ^ G3  =  X2  —  An, 

e^Gi-e^G2  =  (X2-Xi)(G^  +  l), 


Zee, 


(38) 


Proposition  2  The  boundary-value  problem  (|38|)  has  a  unique  solution. 


Proof.  The  proposition  is  equivalent  to  the  fact  that  the  homogenous  problem  (|38[)  has  only  zero  solution, 
which  with  the  representations  ([T4])-([T8])  corresponds  to  (fl3[)  with  the  zero  boundary  conditions,  i.e.  u|5  =  0. 

Let  D+  and  B  be  the  regions  bounded  by  the  surface  S  and  exterior  to  S,  respectively.  The  first  equation 
in  (fT3|)  can  be  recast  in  the  form 

curl£^  +  grad£?+  R^*  —  M2  [k  x  [k  x  (u  -  h— )]]  =0.  (39) 

oz 

With  the  identities 

u  ■  curlu?  =  div(w  X  u)  +  |u>|  , 
u  -  grad  p  =  div(pu)  —  ^?divu  =  div(^?u) , 
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the  scalar  product  of  (39  >  and  u  takes  the  form 


div  (  [wxujf  pn  +  —  |u|2k  )  +  |a;|“  +  M“m,-  ( ur  —  hr  )  =  0. 


On  the  other  hand,  with  the  identity 


[kxh_]-curlh_  =div  ^-|h“|2-h“  (k-h“)  )  , 
the  scalar  product  of  the  second  equation  in  (p~3|)  and  [k  x  h_]  results  in 


K  {K  -«r)--^-div^|h  |2  —  h  (k  h  ))  =0, 


V2 


(40) 


(41) 


(42) 


provided  that  Rm  /  0.  The  linear  combination  (|40])  +  M2([42|)  reduces  to 

div  |[wxu]  +  pn  +  ^|u|2k-  ^^|h^|2  —  h  (k-h  +  |u;|2  +  M2  ( ur  —  h~ )2  =  0. 

Since  curlh+  =  0  in  B+,  we  can  also  write 


(43) 


[k  x  h+l  •  curlh+  =  div 


i  +  |2  _  1,+ 


-h+  (k-h+)  =0. 


(44) 


Let  Hr  be  the  region  bounded  by  the  body’s  surface  S  and  by  a  sphere  Sr  with  large  radius  R  and  center  at 
the  origin,  so  that  D r  — *•  B  as  R  — >  °o.  Applying  the  divergence  theorem  to  the  linear  combination  of  volume 
integrals:  (43  > dV  —  M2/Rm  fffn+  (|44l) dV,  and  using  the  boundary  conditions  u  =  0  and  h+  =  h  on  S, 

we  obtain 


!F  =  Ir  + 


(Jcj|2  +  M2  (ur-h;)2^jdV, 


(45) 


where 


h  =  IL  n'  (  [^xu]+^u+^|u|2k-|^-  (3I11  I2_h  (k  h  )))dS 


Sr 


with  n  being  the  outward  normal.  Note  that  (|45j)  holds  for  multiply  connected  Hr,  since  u,  u,  p,  and  h+  are  all 
single-valued  functions^] 

Next  we  show  that  Ir  — >  0  as  R  — ►  °o.  Let  Cr  be  the  positively  oriented  cross  section  of  Sr  in  the  rz-half 
plane,  i.e.  Cr  is  a  semicircle,  and  let  d/ds  and  d/dn  be  the  tangential  and  normal  derivatives  for  Cr,  respectively. 
With  the  identities 

dr  dz  dr  dz 

ds  dn  ’  dn  dn  ’ 

and  the  surface  element  dS  =  rdsdty,  where  ds  is  the  differential  of  the  length  of  Cr,  we  have 


(46) 


Ir  —  27t  Re 


UcR 


R  \  |\/|2 

(nz  +  iwr)  (  p+iCQ+-(uz-iur)j  +  +  \h;Y  )  rdC, 


(47) 


Now  with  the  representations  (]T4|)-(fT5j)  and  (|TT|)— (p~8|),  the  integral  ([47])  reduces  to 


Ir  —  27t  Re 


UcR 


k]Yx'zG]  +  k2&1KlZG{- 


2A-2  z  2  _ 


(g3“)2 


2Rm(RmR -  M2) 


rdi ; 


5If  Dr  is  multiply  connected,  we  can  make  crosscuts  in  Dr  to  make  Dr  simply  connected,  and  since  u,  w,  p,  and  h±  are  single¬ 
valued,  they  have  the  same  values  on  the  banks  of  a  crosscut. 
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where  k\  =  (R/2  —  2A.2  +  2xA2)  /(k\  —  'k2)2  and  k2  =  (R/2  —  2X\  +  2xAq)  /(Ai  —  A2)2.  Observe  that  Ai  /  A2 
since  M  /  0. 

The  representation  ({8])  implies  that  in  the  spherical  coordinates  (R.  th cp)  related  to  the  cylindrical  coordi¬ 
nates  in  the  ordinary  way,  (7,(7?, if)  =  ,f(tt)R  1  e  A'  R +0  (R  2c  'k:Rj  as  R  — >  00  for  /  =  1,2,  where  _/} (IT J  is  a 
bounded  complex- valued  function  with  |/,('d)|  <  K  for  some  constant  K  and  if  G  [0,7t],  whereas  Example  5  in 
implies  that  G 3  =  0(R~2)  as  R  — >  00.  Thus,  with  ^  =  Re1®,  0  E  [0,7t],  we  can  evaluate 


\IR\  <  2nK2  j  |^.|e-2R(l^l-^C0^)^  sinifdif  +  o(l) 

■4|Xil^  +o(l)  -»0  as  /?-*<», 


y/=i 


where  A,,-  /  0  provided  that  RmR  /  M2.  Consequently,  since  ‘J  =  0,  passing  i?  to  infinity  in  (45),  we  obtain 
///»-  (m2  +  M2(nr-/i7)2)  7F  =  0  so  that  to  =  0  and  nr  =  /ir  in  D  ,  which  along  with  the  representations 

fi~4|),  ([15]),  and  (fl8|)  imply  that  I m Gj  =  0,  j  =  1,2,  in  T)  .  It  follows  from  the  system  ([3])  that  Re Gj  =  Cje~^Jz, 
i  =  1,2,  where  a  are  real- valued  constants.  But  since  G 1  and  G2  vanish  at  infinity,  cj  =  0,  j  =  1,2,  and  thus, 
G\  =  0  and  G2  =  0  in  D  .  Then,  the  first  and  second  equations  of  the  homogeneous  problem  ([38])  imply  Gj  =  0 
and  G3  =  0  on  t,  and  it  follows  from  the  Cauchy  integral  formula  ([5])  for  r-analytic  functions  that  Gf  =  0  in 
©  ,  respectively.  □ 

The  next  proposition  is  considered  from  the  mathematical  point  of  view  only.  It  will  be  used  in  determining 
homogeneous  solutions  to  integral  equations  that  follow  from  the  boundary-value  problem  (|38]). 

Proposition  3  (homogeneous  conjugate  boundary- value  problem,  case  (a))  Let  G 1  and  G2  be  H -analytic 
functions  in  and  let  Gf  be  r-analytic  functions  in  ‘L>+ ,  respectively,  with  G  vanishing  at  infinity.  Under 
the  assumptions  M  f  0,  Rm  f  0,  and  Rm  R  f  M2,  the  homogeneous  conjugate  boundary-value  problem 


1  -  2A2x)  eXl  1 G 1  -  ( 1  -  2Ai  x)  eA2 1 G2  - 


Z, 


—  X2 

RmR -  M2 


G+=0, 


=3iz. 


G\  —  e^2'G2  —  (A.2  —  Ai)G3  =  0, 


^7, 


has  the  solution 


Gl=ce-XiZ,  G2  =  c  e~A2Z,  Gf=2Rmc,  Gf  =  0, 
where  c  is  an  arbitrary  real-valued  constant,  and  Ai  and  X2  are  defined  in  Theorem  [2] 


(48) 


(49) 


Proof.  The  formulas  ([T4])— (fT8]),  in  which  G3  and  G3  are  interchanged  and  G 1,  G2,  and  Gf  are  determined 
in  this  proposition,  represent  u,  fp,  u>,  and  h+  that  satisfy  ([13])  in  D+  and  represent  h  that  satisfy  CD  in 
O-.  In  this  case,  ([13])  in  D+  and  (jTT|)  in  will  be  called  conjugate  equations  for  original  (fT3]l  and  (jTT|). 
The  homogeneous  boundary-value  problem  ([48])  corresponds  to  the  conjugate  equations  with  the  boundary 
conditions  u  =  0  and  h+  =  h  on  S,  where  h  vanishes  at  infinity.  For  the  conjugate  equations,  the  proof  of 
Proposition [4] can  be  repeated  so  that 


F  =  Ir  + 


(Jcj|2  +  M2  {ur-hf)2^dV  =  0, 


(50) 


where 

lR=-tfin{^f-^{k'hb)ds- 
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As  in  the  proof  of  Proposition [4]  it  can  be  shown  that  \Ig\  — >  as  R  — >  °o,  and  passing  R  to  infinity  in  <[50|),  we 
have  u)  =  0  and  ur  =  hf  in  D+,  and  the  representations  (p~4|)— (p~8|)  with  interchanged  G^  and  Gf  imply  that 
I  m  G  i  =  0  and  I  in  Cn  =  0  in  CD+,  and  it  follows  from  the  system  ([3])  that  Re  Gj  =  eye  '~jz,  i  =  1,2,  where  c; 
are  real- valued  constants.  Then  the  second  equation  in  (|48]>  implies  that  G^  =  (<q  —  C2)  / (X2  —  X\ )  on  £,  and 
it  follows  from  the  Cauchy  integral  formula  Q  for  r-analytic  functions  that  Gj  =  (ci  —  cf) / (X2  —  Xi)  in  CD~, 
and  consequently,  G 3  vanishes  at  infinity  only  if  cj  =  C2,  whence  G-,  =  0  in  D  .  Finally,  the  first  equation  in 
([48])  implies  that  G3  =  2Rmc  on  L  where  c  =  c\  =  C2,  and  by  the  Cauchy  integral  formula  ([5])  for  r-analytic 
functions,  we  obtain  G 3"  =  2Rm  c  in  CD+.  □ 

The  problem  (|38j)  can  be  reduced  to  integral  equations  for  the  boundary  values  of  G 1  and  G2  based  on  the 
Cauchy  integral  formula  for  generalized  analytic  functions. 


Theorem  4  (integral  equations,  case  (a))  Let  M  f  0,  Rm  /  0,  and  RmR  /  M2.  In  this  case,  ([38])  yields  two 
integral  equations  for  the  boundary  values  of  Fk(Q  =  e^kZGk(Q,  k  =  1,2: 


(51) 

where  C,  £  (:  and  X]  and  Xj  are  defined  in  Theorem  [2]  A  solution  to  ©  is  determined  up  to  a  real-valued 
constant  c,  i.e.  F/,(Q  =  c,  k  =  1,2,  are  a  homogeneous  solution  to  <HD-  Let  Fk(Q,  k  =  1,2,  solve  then 
Gi(Q  and  G2(Q  on  £  are  determined  by 

Gk{Q  =  (Fk{Q-c)e~**z,  k  =  1,2,  (52) 


where 

c  =  lFk(Q  +  ^~4  Fk( x)e^-Zl)  WH(t;,z,Xk)dx,  C  €  L 

2  27cv/^y^ 


(53) 


Proof.  The  necessary  and  sufficient  conditions  for  G\  (Q,  GjiQ,  and  G±(Q,  £,  to  be  boundary  values  for 

the  corresponding  generalized  analytic  functions  are  given  by  the  generalized  Sokhotski-Plemelj  formulae  ([6]): 


]-Gk(Q  +  ^(j)  Gk{x)WH&,x,Xk)dx  =  0,  F  =  1,2,  (54a) 

2  27ti  Jdjf, 

\g±(Qf^-:4  Gf(z)Wr^,z)dx  =  0.  (54b) 

Z  Z7U  Joj£i 


Expressing  G 3  (£)  and  G-,  (Q  from  ([38])  in  terms  of  G\(Q  and  G 2(Q.  then  substituting  them  into  corre¬ 
sponding  equations  in  (|54b[)  and  solving  the  latter  for  the  integral  terms,  we  have 


Then  (|54a[)  is  multiplied  by  e^K  and  is  recast  in  terms  of  Fk(Q: 


+  4  Fk(z)ek^-^‘WH&,z,Xk)dz  =  0,  F  =  1,2. 

2  2711  ft  U£f 

Subtracting  this  equation  from  the  previous  one  for  corresponding  k,  we  obtain  ([51]). 

Now  we  will  show  that  the  solution  to  the  boundary-value  problem  (|38|)  is  determined  by  ([52])  and  ([53]) 
and  that  a  homogenous  solution  to  ([5l])  is  Fk(Q  =  c,  k  =  1,2.  It  is  known  that  if  a  boundary-value  problem 
for  ordinary  analytic  functions  is  reduced  to  integral  equations,  then  a  homogeneous  solution  to  those  integral 
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equations  is  a  solution  to  the  homogeneous  conjugate  boundary- value  problem;  see  llT8l.  Thus,  the  approach  is 
to  reduce  ([37])  to  the  homogeneous  conjugate  boundary-value  problem  (|48]>. 

Let  Fk(Q,  k=  1,2,  solve  ([37]).  Then  let  0^(Q,  k  =  1,2,  be  //-analytic  functions  in  D+  determined  by  the 
generalized  Cauchy-type  integral 

®k(Q  =  4^4  Fk(x)e-^WH(^h)dx,  *=1,2,  (55) 

^  Jeui' 

and  let  ^(Q  be  r-analytic  functions  in  T)+ ,  respectively,  also  determined  by  the  generalized  Cauchy-type 
integrals 

<L+(Q  =  j)  ((1  -2A,1x)F2(t)-(1-2X2x)F1(t)+X2  — Xj)  Wr(^x)dx,  ?  G  ®+, 

<*>"(«  =  2^-^  (f2(t)  — Fi(x))  ^r(?,T)dT,  C  €  2T . 

With  the  introduced  functions  0+,  02  ,  and  4>  and  the  corresponding  generalized  Sokhotski-Plemelj  formulas 
([6])  for  these  functions,  the  difference  of  equations  ([37])  for  k  =  I  and  *  =  2  is  given  by 

eX2Z0+(Q-ex^0+(Q=d>-(Q, 

Similarly,  the  linear  combination  (1  —  2/i2x)  ■  ©|fc=1-(l-2X1x)-©|  k_2  reduces  to 
(1  -  2Xix)eX2Z0+(Q  -  (1  -2X2x)eXlz0+(Q  =  <f>+(Q,  £  e  t 


Observe  that  the  last  two  equations  are  equivalent  to  the  homogeneous  conjugate  boundary- value  problem 
with  0^  =  Gk,  k  =  1,2,  T>+  =  —  and  d>  =  (Xi  —  A,2)GJ,  Consequently,  by  Proposition  JiJ  the 

only  solution  these  equations  have  is  0^  =  ce~XkZ,  k=  1,2,  <F+  =  2(Xi  —X\)xc,  and  <F  =  0,  where  c  is  a 


real- valued  constant.  In  this  case,  the  representation  (55 )  can  be  rearranged  in  the  form 


1 


4  (Fk(i)-c)e-XkZlWH(^x,h)dT  =  0 

Jeue  v  7 


*=1,2, 


C  G 


2tti  Jfyjf 

For  C,  approaching  l  within  ‘IF  ,  the  above  equation  reduces  to  the  generalized  Sokhotski-Plemelj  formula 


\  (f*(Q  -c)  e~XkZ+^<f^  (f*(t)  -c)  e~XkZl  WH(^T,h)dx  =  0  k=  1,2,  ?  G  F 


(56) 


which  is  the  necessary  and  sufficient  condition  for  (^Fk(Q  —  c j  e  XkZ,  G  £,  to  be  the  boundary  value  of  an 

//-analytic  function  in  ‘D  that  vanishes  at  infinity.  Thus,  the  solution  ([32])-(|33])  follows  from  ([56]). 

Finally,  let  Fk(Q,  k  =  1,2,  £  €  £  be  another  solution  to  ([37]).  Similarly,  we  can  show  that  Gk(Q  = 
£  G  £,  solve  the  boundary-value  problem  ([38]),  where  c  is  a  real-valued  constant. 

However,  by  Proposition [2J  ([38])  has  a  unique  solution,  and  consequently,  Fk(Q  —  Fk(Q  =  c  —  c,  k  =  1 , 2,  ^  €  £, 
which  means  that  a  solution  to  the  integral  equations  ([37])  is  determined  up  to  a  real-valued  constant,  and  the 
proof  is  finished.  □ 

Several  remarks  are  in  order. 


Remark  1  (logarithmic  singularity)  The  kernels  in  ©  have  logarithmic  singularity. 
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Remark  2  (conic  endpoint)  If  i  is  piece-wise  smooth  ( has  salient  points),  in  particular  has  conic  endpoints, 
the  integral  equations  ©  hold  for  all  ^  E  I  except  for  salient  points  and  conic  endpoints,  i.e.  almost  every¬ 
where. 

Remark  3  (multi-connected  body)  The  integral  equations  ©  hold  for  multiply  connected  regions,  e.g.  in 
the  MHD  problem  for  a  torus. 

The  integral  equations  (|5T|)  can  be  solved  by  the  quadratic  error  minimization  method.  Let  i  be  parametrized 
by  £  =  l l(t)  E  C1  [ —  1 , 1],  and  let  Fk  =  Fk(t),  t  E  [0, 1],  k  =  1,2  be  unknown  boundary  values.  For  brevity,  the 
integral  equations  (f5T|)  are  represented  by 

MFl,F2)=fk,  t  E  [—1,1],  k=  1,2, 

where  Ff  arc  corresponding  integral  operators,  and  f  \  it)  =  f2(t)  =  —  l/(2x).  The  functions  Fk(t)  can  be 
approximated  by  finite  functional  series 

m  =  I2^  +  ibk)Tk-i(t),  F2{t)=Y2=lCkTk{t)  +  idkTk_l{t),  (57) 

where  Tk(t)  is  the  Chebyshev  polynomial  of  the  first  kind  and  ak,  bk,  ck,  and  dk,  k  =  1, . . .  ,n,  are  real-valued 
coefficients.  Observe  that  the  real  part  in  the  series  approximating  F2  contains  no  constant  term,  whereas  the 
series  for  F2  does.  This  is  because  same  real-valued  constant  in  places  of  F\  and  F2  is  a  homogeneous  solution 

to  <|5T>. 

Unknown  coefficients  ak,  bk,  ck,  and  dk,  k  =  1 can  be  found  by  minimizing  the  quadratic  error  in 
satisfying  (|5T|>,  i.e. 

min  y“=i||^.(F1,F2)-/,||2  (58) 

ak,bk,ck,dk 

with  the  inner  product  and  norm  for  complex-valued  functions  f(t)  and  g(t)  introduced  b\{^] 

(f,g)  =  R^J  J(t)~^)dt\,  ||/||  =  yftfj).  (59) 

Since  FL/C,  k  =  1,2,  are  linear  operators,  the  problem  (|58|)  is  unconstrained  quadratic  optimization  and  reduces 
to  a  system  of  linear  algebraic  equations  for  finding  ak,  bk,  ck,  and  dk,k=  I , . . .  ,n. 


3.2  Case  (b):  M  ^  0  and  Rm  =  0 

Corollary  [T] implies  that  for  M  f  0  and  Rm  =  0,  the  boundary  conditions  (jTOj)  reduce  to  the  problem  for  deter¬ 
mining  boundary  values  of  G\  and  G2 


ex,zG\  ~eX2ZG2  =  X2~h:  ?  G  £, 

where  both  G\  and  G2  vanish  at  infinity. 


(60) 


Proposition  4  The  boundary-value  problem  ( 60 1  has  a  unique  solution. 


Proof.  The  proof  is  similar  to  the  proof  of  Proposition [4]  In  this  case,  Corollary  [T]  shows  that  the  disturbance 
of  the  magnetic  field  in  and  out  the  body  is  zero,  i.e.  h+  =  0.  With  h  =  0,  the  functional  ([45])  and  the  surface 
integral  (|47|)  reduce  to 


7 


Ir  + 


Dr 


u>\2  +  M2uz 


dV, 


(61) 


6These  functions  are  viewed  as  two-dimensional  vector-functions  from  the  direct  sum  X2([—  1, 1])  ©X2([—  1, 1]). 
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and 


Ir  —  271  Re 


UcR 


(uz  +  iur)  [  fc>+i(Q+^(uz  —  iur)  )  )rdl> 


respecdvely.  Using  the  representation  ([29]),  we  obtain 

r  U2X'zG]-Q1X2ZGl)  cl^ 


2n 

Ir  =  » - —  Re 


Xl  —  X2 


UcR 


and  the  rest  of  the  proof  is  analogous  to  the  proof  of  Proposition  [4] 


□ 


Proposition  5  (homogeneous  conjugate  boundary- value  problem,  case  (b))  Let  G\  and  G2  be  H -analytic 
functions  in  'D+ .  Under  the  assumption  M  f  0,  the  homogeneous  conjugate  boundary-value  problem 


:Gi-e^G2  =  0,  £  G 


has  the  solution 


G  i  =ce4lZ, 


G2  =  ce4jz, 


(62) 

(63) 


where  c  is  an  arbitrary  real-valued  constant,  and  X]  and  X2  are  defined  in  Corollary  [7] 


Proof.  The  proof  is  similar  to  the  proof  of  Proposition  [3]  □ 

With  the  generalized  Cauchy  integral  formula  for  //-analytic  functions,  the  boundary-value  problem  (T0|) 
readily  reduces  to  an  integral  equation. 

Theorem  5  (integral  equations,  case  (b))  Let  M  f  0  and  Rm  =  0,  and  let  F]  (Q  =  e^lzGi( Q.  The  boundary- 
value  problem  d60])  reduces  to  the  integral  equation  for  the  boundary  value  ofF\: 


1  i  Lxdz-zd  xM))  Ffx)dx  =  X2-XU  C  G  l 

27t  t  Je\j£i  \  / 


(64) 


A  solution  to  ( 64 1  is  determined  up  to  a  real-valued  constant  c.  Let  F\  (Q  solve  ( 64 1,  then  G\  (Q  in  ( 60 1  is  given 
by 

G1(Q  =  (fi(Q-c)e-^z,  ?€*,  (65) 


where 


c  =  +  4  T1(x)e^(^)  WH&,x,Xi)dx, 

2  2m 


?€• 


(66) 


Proof.  The  integral  equation  ([64])  is  derived  similarly  to  ([57]).  In  the  second  part  of  the  proof  that  determines 
a  homogeneous  solution  to  (|64j),  //-analytic  functions  0j| ,  k  =  1,2,  in  D 1  are  introduced  by  the  generalized 
Cauchy-type  integrals 

=  Fi(x)e-liZ1WH^,x,Xk)dx,  U®+, 

2rct  Jtue> 


&2(Q  =  ^<fe  ii{Fl(x)+Xl-X2)^-X2Zl^H^,xM)dx,  U®4 


where  F\  is  a  solution  to  (64 1.  With  these  functions  and  corresponding  generalized  Sokhotski-Plemelj  formu¬ 
las  ([6]),  the  integral  equation  ([64])  reduces  to  the  homogenous  conjugate  boundary-value  problem  (|62|),  whose 
solution  is  given  by  ([63]).  The  rest  of  the  proof  is  analogous  to  that  of  Theorem[4]  □ 

The  integral  equation  ([64])  can  be  solved  by  the  quadratic  error  minimization  method  outlined  in  the  end  of 
Section|3T|  Remarks [T|-[3] also  hold  for  d64j). 
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3.3  Case  (c):  M  ^  0  and  RmR  =  M2 

For  M  /  0  and  RmR  =  M2,  the  boundary  conditions  (fTO])  reduce  to  the  boundary-value  problem  for  the  four 
generalized  analytic  functions  Gi,  G2,  and  G±  defined  in  Theorem  [3] 

2e^G1-G2  +  (R  +  Rm)(G+  +  l)=0,  £  e 
Rm((z-^-)G2  +  G3-  +  l) -RG+  =  °,  t()'! 

where  A  =  (R  +  Rm)/2. 

Proposition  6  The  boundary -value  problem  (|67[)  has  a  unique  solution. 


Proof.  The  assumptions  M  /  0  and  RmR  =  M2  imply  Rm  /  0,  and  consequently,  as  in  the  proof  of  Proposition 
[4j  we  obtain  the  functional  (|45|)  with  the  surface  integral  (|47|).  With  the  representations  ([30])— ([34]),  the  integral 
<@7)  reduces  to 


Ir  =  2n  Re 


uz  -j-iur)G2  + 


R 

2(R  +  Rm) 


(2eXzGi 


Since  in  this  case,  G\  and  G2  are  //-analytic  and  r-analytic  functions,  respectively,  we  have  G\  (R.i})  = 
/(fi)/?_1  e~^R  +0  as  R  — >  00,  where  /(fi)  is  a  bounded  complex-valued  function,  and  G2  = 

0(/C2)  as  R  —>  00.  The  rest  of  the  proof  is  completely  analogous  to  the  proof  of  Proposition |dj  □ 


Proposition  7  (homogeneous  conjugate  boundary-value  problem,  case  (c))  Let  G 1  and  G2  be  H -analytic 
and  r-analytic  functions  in  CD+,  respectively,  and  let  G±  be  r-analytic  functions  in  respectively,  with 
G  vanishing  at  infinity.  Under  the  assumptions  M/0  and  RmR  =  M2,  the  homogeneous  conjugate  boundary- 
value  problem 

2e^Gi-G2  +  (R  +  Rm)G3-  =  0,  Z  e  l, 

Rm((z4r)G2  +  G+)-RG3=°,  tei,  (68) 

has  the  solution 

Gi  =  ce-Xz,  G2  =  2c,  G+  =  —  c(2z-ir),  Gf  =  0,  (69) 

where  c  is  an  arbitrary  real-valued  constant,  and  A  =  (R  +  Rm) /2. 

Proof.  The  proof  is  similar  to  the  proof  of  Proposition  [3]  □ 


Theorem  6  (integral  equations,  case  (c))  Let  M  /  0  and  RmR  =  M2,  and  let  A  =  (R+  Rm)/2,  then  the  boundary - 
value  problem  (|67j)  reduces  to  two  integral  equations  for  /q(Q  =  e^Gi(Q  and  F2(  Q  =  G2(Q,  Z  F  F. 


(e^"Zl)  ~  FiWdt  +  ftiZ)  -  \f2{Z)  =  0, 

(zi -z+ i(r-ri0  mr%t(2FiK) 


R+  Rm 


(70) 


Equations  ©  determine  F\  and  F2  up  to  c  and  2c,  where  c  is  a  real-valued  constant.  Let  /),  k  =  1,2,  solve 
©  then  the  solution  to  the  boundary-value  problem  m  is  given  by 

Gl(Z)=(Fl(Z)-c)e-Xz,  G2(Z)  =  F2(Z)  —2c,  (71) 

where 

2  c=\hQ  +  ^~4  F2(z)Wr(Z,z)dz,  ZfL  (72) 

2  27U  fe\jp 
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Proof.  The  derivation  of  (|70|)  is  similar  to  that  of  (|57|);  see  the  proof  of  Theorem [4]  In  the  second  part  of  the 
proof  that  determines  a  homogeneous  solution,  an  //-analytic  function  0+  in  ‘D ",  r-analytic  functions  <T> ,  and 
(P2  in  ‘D  ,  and  an  r-analytic  function  <P2  in  (D~  are  introduced  by  the  generalized  Cauchy-type  integrals 

®+(Q  =  ^4 

27tt  j£ \jp 

®t(Q  =  4fi4  F2(%)Wr^,x)dX,  ?€©+ 

((21  -  u) *w+ mrW  (*«  -sm) + ^=)  ®+> 

^2  (0  =  2^-^  (2Ti  (x)  -F2(t))  Wr&,x)d?,  ?  e  2T , 

where  Ft,  k  =  1 , 2,  are  a  solution  to  (|70|).  With  these  functions  and  corresponding  generalized  Sokhotski-Plemelj 
formulas  ([6]>,  the  integral  equations  (|7()|  reduce  to  the  homogeneous  conjugate  boundary-value  problem  ([68]). 
whose  solution  is  given  by  (|69|).  The  rest  of  the  proof  continues  as  in  the  proof  of  Theorem[4j  □ 

The  integral  equation  (fTO])  can  be  solved  by  the  quadratic  error  minimization  method  outlined  in  the  end  of 
Section |3T[  Remarks [T]-[3] hold  for  (f70j)  as  well. 


3.4  Drag 

The  drag  (force)  exerted  on  the  body  by  the  electrically  conducting  fluid  in  the  presence  of  the  magnetic  field 
has  the  mechanic  and  magnetic  components 

F  =  ff  P ndS-n  [If  [J  x  H]  dV, 

JJ  s  JJJd- 

where  P„  =  pv  (2 du/d/;  +  [n  x  curl u] )  —  pn,  with  n  being  the  outward  normal  to  S,  and  J  =  d-  curl H. 

In  the  axially  symmetric  case,  divu  =  0  and  the  boundary  conditions  ([TO])  imply 

=  — [nxcurlul  on  S,  (73) 

an 

see  e.g.  the  proof  of  Proposition  8  in  ll28l.  Also,  the  second  term  in  F  vanishes. 

Proposition  8  If  at  infinity,  the  flow  and  magnetic  field  are  constant  and  parallel  to  the  body’s  axis  of  revo¬ 
lution,  and  if  also  the  fluid  and  body  have  the  same  magnetic  permeability,  the  magnetic  force  has  no  direct 
contribution  to  the  drag. 


Proof.  In  the  axially  symmetric  case,  both  the  total  mechanic  and  magnetic  forces  are  parallel  to  the  body’s 
axis  of  revolution,  and  consequently,  it  is  sufficient  to  show  that  the  projection  of  the  magnetic  force  onto  k 
vanishes.  Let  D r  be  the  region  bounded  by  the  body’s  surface  S  and  sphere  Sr  with  large  radius  R  and  center  at 
the  origin.  Using  the  identity  (|4T|)  and  the  divergence  theorem,  we  obtain 


k  •  [J  x  H]  dV 


-^JJJJ‘x'rhcaAh~dv 


jj 2 

7) - 


where 


k  h~)h  —  ^|h  |2k  JdS, 


IR=  /I  n-  (  (k-h_)h“  —  ^|h“|zk  )  dS. 

Sr 
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As  in  the  proof  of  Proposition [4]  we  can  show  that  in  all  three  cases  (a),  (b),  and  (c),  Ir  — >  0  as  R  — >  °o.  On  the 
other  hand,  since  h+  =  h  on  S  and  curlh+  =  0  in  B+,  using  the  identity  similar  to  (|4T|)  for  h+,  we  have 


/=  I/  n- (  (k-h+)h+--|h+|2k)rfS=- 


[k  x  h+l  •  curl  h+  dV  =  0. 


Consequently,  fffD  k  •  [J  x  H]  dV  — >  0  as  R  — >  °o.  □ 

In  view  of  (f73|)  and  Proposition [8j  the  drag  in  terms  of  dimensionalized  u  and  p  takes  the  form 

(74) 


Fz  =  —  Voopva  j  j  k  •  ([n  x  curlu]  +  pn  )dS, 


which,  as  in  the  proof  of  Proposition  1 1  in  If28l ,  reduces  to 


Fz  =  -27tVk,pvaRe  r(p+ioi)d^ 


(75) 


Let  Cd  =  —Fz/(6nVoopva)  be  a  dimensionless  drag  coefficient,  where  67t\4=pv«  is  the  drag  of  the  sphere 
with  radius  a  in  the  Stokes  flow.  Theorem[2]  Corollary[l]  and  Theorem[3]and  the  corresponding  boundary-value 
problems  (|38|),  (|60|),  and  (|67|)  imply  that 


p  +  ico=  < 


1 


Xj  —  Xi 

Ai  z, 


(R-2^2)eXlZGi-(R-2^1)e5l2ZG2j  +  R,  case  (a), 

2(e^Gi+Xi), 

1 


2ReXzGi  +  RmG2)  +  R, 


C,  e  £  case  (b), 
£  G  i  case  (c). 


,  R  +  Rn 

Now  with  this  representation  for  p+ ico  and  the  fact  that  Re  [J(rdQ  =  0,  ([75])  leads  to  the  following  result. 

Proposition  9  (drag)  Let  M  /0,  then  the  drag  coefficient  Cd  is  determined  as  follows. 

(a)  If  Rm  7^  0  and  RmR  p  M2,  then 

1 


Cd  = 


Re 


3(Xi  —  Xf) 

where  G\  and  G2  are  defined  in  Theorem  [2] 
(b)  If  Rm  =  0,  then 


J^r  ((R  -  2X2)eXlZGi  -  (R  -  2Xi)eX2ZG2)  </£ 


(76) 


CD  =  3  Re 


[  re^lZGi  dC, 

Jt 


(77) 


where  G i  and  K\  are  defined  in  Corollary^ 7] 
(c)  7/RmR  =  M2,  then 


CD  = 


3(R+R 

where  G\  and  G2  are  defined  in  Theorem  [J] 


Re 


Jr  ^2R  e^zGi  +  RmG2^  d £ 


(78) 


19 


4  Axially  Symmetric  MHD  Problems 

4.1  Sphere 

The  purpose  of  solving  the  MHD  problem  for  sphere  is  three-fold:  (i)  verifying  solution  representations  in 
terms  of  generalized  analytic  functions  by  comparing  to  the  existing  solutions  for  sphere;  (ii)  testing  accuracy 
of  solutions  to  the  integral  equations,  and  (iii)  obtaining  drag  for  various  R,  Rm,  and  M  to  compare  to  the  drag 
of  minimum-drag  shapes. 

Let  Sa  be  the  sphere  with  radius  a  and  center  at  the  origin,  and  let  (A'.thcp)  be  the  spherical  coordinates 
related  to  the  cylindrical  coordinates  in  the  ordinary  way. 


4.1.1  Case  (a):  M  ^  0,  Rm  ^  0,  and  Rm  R  /  M2 

In  the  region  exterior  to  Sa,  the  functions  G\  and  Gi  in  Theorem[2]can  be  represented  by  series  similar  to  ([8]): 


Gk(R,ft)  =  \  ^  Y,Ak,nNn(cos$,R,Xk),  k  =  1,2, 


(79) 


n—\ 


where  Ak  n  are  unknown  real-valued  coefficients  and 

k„+.(i  m 


N„(t,R,l)=Ln(t) 


-  (signX)  L_„(f) 


KJ(_i  (|A,|i?) 

K„+i(Wa)’ 


whereas  Gj  can  be  represented  by 


Gf(R,*)  =^=^{R/a) ^L^cosfl), 


(80) 


where  B±  are  unknown  real- valued  coefficients;  see  lf30l  Example  5]. 

With  the  orthogonality  property 

(L,;(f)  ,Lm(f))  =  2/7i  8mn  (81) 

for  all  integer  n  and  m,  where  8mn  is  the  Kronecker  delta,  and  the  representations  (|79|)  and  (80 ».  the  inner  product 
of  the  first  equation  in  (f38|)  with  L_m_i(f)  for  m  >  0  and  the  inner  product  of  the  second  equation  in  (|38j)  with 
L  m(t)  for  m  >  1  reduce  to  a  real- valued  infinite  linear  system  for  A\n  and  Ajy. 

(1  —  2^2x)J^n=l  ~  (1  ~  kmn{^2)M,n  =  2(^1  ~  ^2)§m0>  m  >  0, 


)-Al,n  ^innC^-l)  ^2,n  — 0, 


m  >  1, 


(82) 


where  =  (e^at  N„(f,aA),L_m_i(f))  and  =  (eXat  Nn(t,a,X),Lm{t)).  The  system  ([82])  can  be 

truncated  at  some  large  m  and  solved  numerically. 

With  the  second  equation  in  ([82]),  the  drag  coefficinet  ([76])  reduces  to 

2  oo 

CD  =  ^  X2n=l^1"rl1”(^“l)- 

Figure  [T]  shows  Cd  as  a  function  of  S  for  five  pairs:  R  =  Rm  =  1;  R  =  1,  Rm  =  3;  R  =  Rm  =  2;  R  =  3, 
Rm  =  1;  and  R  =  Rm  =  3. 
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4.1.2  Case  (b):  M  /  0  and  Rm  =  0 


As  in  the  previous  case,  the  functions  G i  and  (A  in  Corollary  [T] can  be  represented  by  the  series  f79]).  With  the 
orthogonality  property  (f8T[)  and  representations  ([79]),  the  inner  product  of  ([60])  with  e  'k2alLm(t)  for  m  >  1  and 
the  inner  product  of  (|60[)  with  c  kiaILm(t)  for  m  >  1,  reduce  to  a  real-valued  infinite  linear  system  for  A\n  and 

A2  ,«• 

|  '^mn  A  1 .11  2  HI /\ 2jn  ~  ^2. in  -  777  7  1  , 

^-,00  „  (w) 

2l7lA\m  _ |  ffiwz ^ 2,n  —  b\ Wl  A  1, 


where  in  this  case,  t,„ 


=  (eXiafN„(t,a,^i 


Lm(t)),  T]m„  =  (eX2a?N„(f,a,X2),e  XlflfLm(f)),  and  bk,m  = 


(X2  —  'hi,e~^katl-m(t)),  k  =  1,2.  Similarly,  the  system  (83)  can  be  truncated  at  some  large  m  and  solved 
numerically. 

The  drag  coefficient  (f77|)  reduces  to 


CD 


(eWN  n{t,aM)M{t))- 


Figure  [2] depicts  Co  as  a  function  of  R  for  Rm  =  0  and  M  =  0, 1,2,3,  and  as  a  function  of  M  for  Rm  =  0 
and  R  =  0, 1,2,3.  Numerical  results  show  that  for  R  =  Rm  =  0,  the  drag  is  in  good  agreement  with  Chester’s 
formula  |9|:  Co  =  1  +  |M  +  g^M2  —  ^^IVF  +  0(M4)  for  small  M  up  to  0.25;  see  curve  a  on  Figure[2j}. 


4.1.3  Case  (c):  M  ±  0  and  Rm  R  =  M2  (S  =  1) 

In  this  case,  the  //-analytic  function  G  i  and  the  r-analytic  function  Gi  in  Theorem  [3]  are  represented  by  the 
series 


/oo 

I  52^t,n N„ (cos'd,/?, A,),  G2(/?,i3)=£;=1A2!n(/?/a)-n-1L„(cosi3), 

n—  1 


(84) 


4.1.1 


and  X  =  (R  +  Rm) /2,  whereas  the  r-analytic  functions  G^ 

With  the  orthogonality  property  ([8T])  and  representations  ([84])  and  ([80]).  the  inner  products  of  the  first  equa¬ 
tion  in  ©  with  L m(t)  for  m  >  1  and  with  L  _m_  i  (t)  for  m  >  0  reduce  to  a  real- valued  infinite  linear  system 


where  N„ (cos'd, R,X)  is  introduced  in  Section 
are  represented  by 


m 4 2jn  — 0,  m  A  1, 

y. _ 1  ^,mn{X)Aitn  —  (R  +  Rm)(/77  +  l)/?jj+i  =  (R+  Rm)8mo,  777  A  0, 


(85) 


where  x\mn(X)  and  are  those  in  Section  4.1.1 


Similarly,  with  ([84])  and  (80  (  and  the  orthogonality  property  for  the  Legendre  polynomials,  the  inner  prod¬ 
ucts  of  the  second  equation  in  ([67])  with  L m(t)  for  m  >  1  and  with  \  (t  )  for  m  >  0  reduce  to 


m—  1 
2m  —  1 


A2.m_i  + 


3  (m  +  1 ) 
2m  +  1 


M,m+ l)+2Sm  =  0,  777  >1, 


Rm  a  ,[’n+  1^'”  +  2L  A2;m+1  -  R(t77  +  l)B+  =  -Rm  8m0,  m  >  0, 


(86) 


‘  (2/77+  1  )(2/77  +  3 "V"*  '  V  m+1 

It  follows  from  the  system  ([85])  and  the  second  equation  in  (86 )  that 

f  fl,Rm(R"F  Rm )  (t/7  +  2)  -  /a  ,  \  , 

V  R(2m  +  l)(2m  +  3)  W  “ U W J  -  - 


R  +  Rn 
R- 


^m0  5 


777  >  0. 


21 


The  drag  coefficient  (|78])  simplifies  to  the  formula 

2  »  1  oo  . 

CD  =  3  Tn=l^1.”Tlln^)’ 

which  is  used  to  evaluate  Cd  at  S  =  1  in  Figure[l]and  which  for  various  small  R  and  Rm  agrees  with  Gotoh’s  for¬ 
mula  m-  CD  «  1  +  |R-  (^,R2  +  BRRm)  +  (4o  (213R3  +  256R2Rm -704RRm2)  +  ^RRm(R  +  4Rm)). 


(a)  curves  a,  b,  c,  d  correspond  to  M  =  0, 1 , 2, 3  (b)  curves  a,  b,  c,  d  correspond  to  R  =  0, 1 , 2, 3 

Figure  2:  Drag  for  the  unit  sphere  normalized  to  the  sphere’s  Stokes  drag  GnvpVoo  for  Rm  =  0  as  functions  of 
R  and  M. 


4.2  Minimum-drag  Spheroids 

This  section  solves  the  axially  symmetric  MF1D  problem  for  solid  nonmagnetic  spheroids  having  the  volume  of 
a  unit  sphere,  i.e.  4tt / 3,  and  among  those  finds  the  spheroids  that  have  the  smallest  drag  for  given  R,  Rm,  and 
M.  In  the  rz-half  plane  (r  >  0),  the  cross  section  i  of  spheroids  is  parametrized  by  r(t)  =  a  cost,  z(t)  =  cT~  sint, 
t  E  [— 7t/2,7t/2],  with  a  E  (0,1].  Then  for  fixed  R,  Rm,  and  M,  the  spheroid  drag  becomes  a  function  of  a, 
and  the  minimum  of  this  function  along  with  the  value  of  a  at  which  the  minimum  is  attained  is  found  by  the 
modified  bisection  method. 

Let  Cd  and  C*D  be  the  drag  coefficients  for  the  minimum-drag  spheroid  and  unit  sphere,  respectively,  for 
same  R,  Rm,  and  S,  and  let  K  be  the  axes  ratio  of  the  corresponding  minimum-drag  spheroid.  Figure [3] shows  the 
drag  ratio  Cd/Cq  and  K  as  functions  of  S  E  [0,2]  for  the  three  pairs  R  =  Rm  =  1,  R  =  Rm  =  2,  and  R  =  Rm  =  3. 
As  a  function  of  S,  Co /C'D  has  maximum  at  S  =  1  (and  is  nonsmooth  at  S  =  1),  which  suggests  that  drag 
reduction  is  smallest  for  S  =  1  and  is  more  significant  for  S  »  1,  whereas  the  shortest  minimum-drag  spheroids 
appear  to  be  not  at  S  =  1 .  The  found  minimum-drag  spheroids  will  be  used  as  initial  approximations  for  the 
minimum-drag  shapes  (Section[8]). 

5  Optimality  Condition  for  Minimum-drag  Shape 

This  section  derives  necessary  optimality  conditions  for  the  shape  of  the  body  that  has  the  smallest  drag  and 
the  volume  of  a  unit  sphere  in  the  MHD  problem  considered  in  Section  [2] 

Let  body’s  surface  S  divide  the  space  into  the  interior  and  exterior  regions  D+  and  D~,  respectively.  Since 
the  drag  F  is  parallel  to  the  z  axis,  minimizing  the  absolute  value  of  F  is  equivalent  to  minimizing  — k  •  F,  and 
in  view  of  (|74|),  the  shape  optimization  problem  is  formulated  by 

mm's,  E  =  k-  jj  ([n  x  curlu]  +  fc>n)dS  (87) 
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(a)  drag  ratio  Cd/C*d  versus  S  (b)  axes  ratio  versus  S 

Figure  3:  The  drag  ratio  Cd/C*d  for  the  minimum-drag  spheroid  and  corresponding  spheroid’s  axes  ratio  as 
functions  of  S  for  R  =  Rm  =  1  (curve  a);  R  =  Rm  =  2  (curve  b);  and  R  =  Rm  =  3  (curve  c). 


subject  to  (fTT|)  and  (p~3|)  reformulated  as 

curlcurlu  +  grad£?+  R-*1  -  M2  [k  x  [k  x  (u-h-)]]  =0, 
oz 

*  curlh- =  Rm  [(u  — h“)  x  k], 
divu  =  0,  divh_=0. 

with  the  boundary  conditions  (fl2])  and  (flO]). 

To  derive  necessary  optimality  conditions  for  (|871)-([88j),  we  use  Mironov’s  shape  variation  approach  ITT7I. 
Let  r  be  the  radius  vector  describing  the  optimal  shape  S,  and  let  the  shape  variation  Se  be  determined  by 


tE  =  t  +  e/(t)n, 


where  £  is  a  positive  small  number,  n  is  the  outward  normal,  and  /( t)  is  any  continuous  bounded  scalar  function 
such  that  jfsfdS  =  0.  The  condition  on  /  follows  from  the  volume  variation:  if  o+  is  the  region  bounded  by 
the  surface  SE,  then 


dV+e  I  I  fdS+o(e)  =4n/3, 


and  since  fffD+  dV  =  4k/3,  f  should  satisfy  JfsfdS  =  0. 

Let  u,  p,  and  h+  be  the  velocity  disturbance,  pressure,  and  magnetic  field  disturbances  in  and  out  the  body 
for  the  optimal  shape,  and  let  their  variations  be  given  by 


uE  =uT£Ui+o(e),  pe  =  p+epi  +o(e),  =u  +  £hjb  +  o(s), 

where  m,  p\,  and  h satisfy  (|88j). 

The  boundary  condition  (jTOj)  should  hold  for  uE  on  5e: 


uE  =  -k  on  Se, 


which  implies 


uE  =  u(r  +  £/(r)n)  +  £Ui(r  +  £/(r)n)  +  o(£)  =  -k. 
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whence 


and,  consequently, 

or,  equivalently, 


u(r)  +  8(/(r)3^)+u1(r))+0(s)  =  -k 


3u 

/  +ui  =  0  on  S, 
on 


ui  =  —  =  f\n  x  curl ul  on  S. 

dn  1  1 


Similarly,  the  boundary  condition  <[T2|)  should  hold  for  on  Se: 


h+  =  he  on  Sg, 


(89) 


which  implies 

h+(t  +  £  /( r)  n)  +  eh|(t  +  e. /( r)  n)  +  o(e)  =  h~  (t  +  £  /( r)  n)  +  (r  +  £  /( t)  n)  +  o(e) , 

whence 

h+(r)  +  £  (f(x)  3*^r)  +  h+(t))  +  o(E)  =  h-  (r)  +  £  (f(x)  3*^r)  +  hr (t))  +  o(e) 


and,  consequently, 


h+-hr  =  -/ 


3h+  3h 


dn  dn 


on  S. 


The  variation  of  the  functional  £  defined  in  ([87])  can  be  determined  as  follows: 

£e  =  k-  //  ([n  x  curluE]  +  pzn)dS 

JJse 

=  k  //  ( [n  x  curl ue]  +  n) dS + k •  //  ([n  x  curluE]  +  pzn)dS  =  /E  +  /E, 

JJ  s  JJ  sz-s 


=4 


=4 


with 


and 


4  =  k  • 


7E  =  £  +  £k-  Jj  ([n  x  curlui]  +  pi  n)dS+o{e), 


(curl  curl  u  +  grad  p+£  (curl  curl  uj  +  grad  p\ ))  sign/ dV  +  o(e) 


du- 


=  £k-  //  (curlcurlu  +  grad^?)/dS'  +  o(£)  =  — eR  //  fdS +  o(e). 


JJs 


dz 


(90) 


(91) 


(92) 


(93) 


where  DE  is  the  region  bounded  by  the  surfaces  S  and  SE,  and  in  obtaining  (|93|),  we  used  the  divergence  theorem, 
the  first  equation  in  (|88j)  and  the  fact  that  k  •  [k  x  [k  x  (u  —  h~]]  =  0. 

A  further  transformation  of  the  functional  /E  relies  on  the  equations  adjoint  to  i 


Proposition  10  (adjoint  equations)  Let  w  G  C2(D  ),  q  G  C]  (D  ),  g1*1  G  and  v±  G  satisfy  the 

adjoint  equations 

curl  curl  w-  R— - M'[kx[kx  w]]  +  grader  +  Rm  [g_  x  k]  =  0,  divw  =  0, 

dz  (94) 

curlg~  +  gradv“  =  Rm  [g_  x  k]  -  M2  [k  x  [k  x  w]], 
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in  D  and 

curlg+  +  gradv+ =  0  in  IB f 

with  the  boundary  conditions  on  S  and  conditions  at  infinity  to  be 

,+  |  — 


Then 


wL  =  k,  g+L  =  g+|c, 

w|oo  =  0,  =  0,  v“|oo  =  0,  ?U  =  o. 


k-  //  ([n  x  curlui]  +  p\n)dS  =  //  ui  •  ([n  x  curlw]  +qn—  R(n-  k)k)<LS 


+  jjn- ([(h^-hj )  xg  ]-(h+-h!  )v  )dS, 


provided  that 


Ir=  I  n- ([curluj  x  w]  +  [ui  x  curlw] +/?iw  —  gui 

lsR 


+R(w-ui)k  +  [h1  xg  ]— hj  -v  )dS— >0  as  0) 

where  Sr  is  the  sphere  centered  at  the  origin  and  having  large  radius  R. 


Proof.  Since  m,  p\,  and  hj  satisfy  ([88]),  we  can  write 

0  =  w  ^curlcurlui  +  R^-  -  M2  [k  x  [k  x  (m  -  hj-)]]  +  grad pj  —  gdivui 

+  g_  •  (curlhj"  —  Rm  [(uj  —  h]f)  xk])-v"  divhj~ 

=  div  ([curlui  x  w]  +  [m  x  curlw]  +  piw  —  qui  +  R  (w-ui)kd-  [hj"  x  g_]  —  •  v_) 

+  ui  •  ^curl curlw  —  Ry  -  M2  [k  x  [k  x  w]]  d- graded-  Rm  [g~  x  k]^  —p\ divw 

d-hj'  •  (curl  g“  d- grad  v“  -  Rm  [g_  xk]  +  M2[kx[kx  w]]) , 

which  in  view  of  the  adjoint  equations  ([94])  implies 

div  ([curlui  xw]  +  [ui  x  curlw]  d- p\v/  —  gui  -f  R(w-ui)k-f  [hj-  x  g-]  — hj"  •  v“)  =  0  in  ID 

On  the  other  hand,  since  h+  satisfies  curl  h+  =  0  and  div  h+  =  0,  we  have 

0  =  g+ -curl  —v+ div  hj1"  =  div([h+  x  g+]  — v+h^)  d-h^  •  (curl  g+d-  grad  v+) , 


which  with  the  adjoint  equation  ([95])  reduces  to 

div  ([hj"  x  g+]  —  v+h]H)  =  0  in 


(95) 


(96) 


(97) 


(98) 


(99) 


Let  D/}  be  the  region  bounded  by  the  surface  S  and  sphere  Sr.  Integrating  ([98])  and  (|99[)  over  Hr  and  B+, 
respectively,  and  then  using  the  divergence  theorem  and  adding  the  resulting  equations,  we  obtain 


Ir  —  JJ  n-  ([curlui  xw]  +  [ui  x  curlw]  d- piw  —  qu\  +  R(w-ui)k 
+[K  xg"]-h"-v"-[h|  xg+]d-v+h+)dS  =  0, 

which  with  the  boundary  conditions  (|96|)  and  the  assumption  that  Ir  —  0  as  R  — >  °o  can  be  rearranged  in  the 
form  d97[).  □ 
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With  the  relationship  ([97])  and  the  boundary  conditions  (89 )— ( 90 ),  the  functional  (92)  takes  the  form 
Ie  =  ‘E  +  z  Jj  ^[n  x  curl u]  •  [n  x  curl  w]  +  R(n-k)^^  fdS 


~ e'lsn ' 


3h+  3h 


dn  dn 

Now  since  du/ds  =  0  on  S,  we  observe  that 


xg 


dh+  dh 


dn  dn 


(100) 


v  I  fdS  +  o(e). 


du- 

dz  duz 

dz  j 

(  du-  dr  du-  dz  \ 

dr  I 

(  du-  dz 

du-  dfj 

.  du- 

dn 

dn  dn 

dn  ' 

c  dr  dn  dz  dn  J 

dn  ' 

\  dr  dn 

dz  dn  J 

^  +  ~d^ 

dr  duz  duz  duz 
dn  ds  dz  dz 


(101) 


Let  h  =  h+  —  on  S.  Then  (fT2[)  and  the  conditions  divh^  =  0  in  D1*1  U  S  imply  that  h  =  0  and  divh  =  0 
on  S.  In  this  case,  dh/dn  =  —  [n  x  curl/?]  on  S',  see  the  proof  of  Proposition  1 1  in  1281.  On  the  other  hand,  it 
follows  from  (fTT|),  the  second  equation  in  ([13]),  and  (p~0[)  that  curl/?  =  Rm  [h  x  k]  on  S.  Thus, 


dh+  dh 


dn  dn 


=  -Rm[nx[h  xk]]  on  S. 


(102) 


Also,  in  the  axially  symmetric  case,  curlu  =  to (r,z)  e<p,  curlw  =  CO *(r,z)  ejp,  and  g  =g  ( r,z )  e<p. 
Finally,  in  view  of  (fJ3|)  and  (|  1 00|)  with  (jTOT])— (JT02]),  (|9T|)  reduces  to 


(Eg  —  (E  T  £  jj  (co  co  +  Rm  hr  g  ^  f  dS  T  c?(e)  . 

The  necessary  optimality  condition  requires  lim(Ee  —  £)/e  =  0,  which  yields  ffs  (coco*  +  Rm  h~g~)  fdS. 

However,  since  /  should  satisfy  ffsfdS  =  0,  we  conclude  that  coco*  +  Rm  /?,7g-  is  constant  on  S. 
Consequently,  the  following  result  has  been  proved. 

Theorem  7  Under  the  assumptions  of  Proposition  [7d]  the  necessary  optimality  condition  for  the  minimum- 
drag  shape  subject  to  the  volume  constraint  is  given  by 


coco*  +  Rm/?r  g  =  const  on  S, 


(103) 


where  CO  =  em  •  curlu,  CO*  =  em  •  curlw  and  g  =  e(p  •  g  ,  and  w  and  g  are  solutions  of  the  adjoint  equations 


(94 1— (95 )  with  the  boundary  conditions  (96 ). 


6  Solving  Adjoint  MHD  Equations 


A  solution  to  the  adjoint  equations  ([94[)-(|95])  with  the  boundary  conditions  (|96|)  can  be  also  constructed  in  terms 
of  generalized  analytic  functions.  In  the  axially  symmetric  case  with  the  axis  of  revolution  coinciding  with  the 
z  axis,  the  functions  w,  q,  g*,  and  v+  have  the  following  representations: 

w  =  Wr{r,z)er  +  wz(r,z)K  q  =  q(r,z),  g±=g±(r,z)^,  v±=v±(r,z), 

where  the  third  formula  implies  divg+  =  0. 
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6.1  Case  (a):  M^O  and  RmR  ^  M2 

Let  i  be  the  positively  oriented  cross  section  of  S  in  the  right-half  rz- plane,  and  let  (!  be  the  reflection  of  l  over 
the  z  axis. 


Theorem  8  (solution  representation  for  the  adjoint  equations,  case  (a))  Let  X\  and  X2  be  as  those  defined 
in  Theorem  [2]  In  the  axially  symmetric  case  with  M  f  0  and  RmR  f  M2,  a  solution  to  the  adjoint  equations 
494])— (|95|)  is  given  by 


wz  +  iw,.=  ^  ((2X2x-l)e^lzGi-(2Xix-l)e^2ZG2)+xG^, 

v~  +  ig~  =  (1  -  Rx)  (x2e~XlZGl  +  Gj\  , 

v+  +  ig+  =  Gj-,  (104) 

q=2h-x!^  (^-26  ^"uRe[Gi]  —  X,! e~^2ZRe[G2]^  +  Rx Re[Gj], 

C0*  =  X^X2  ((R-2^2)e^i:Im[G1]-(R-2X1)e^22Im[G2])  . 


where  G\  and  G2  are  H -analytic  functions  in  D  that  satisfy  ©  with  X  =  —X\p,  respectively,  and  vanish  at 
infinity;  Gf  are  r- analytic  functions  in  CD±,  respectively,  with  Gj  vanishing  at  infinity. 


Proof.  The  proof  is  similar  to  that  of  Theorem  [2]  □ 

With  the  representations  ( 104[).  the  boundary  conditions  (96 1  reduce  to  the  boundary-value  problem  for  the 
generalized  analytic  functions  G\,  G2,  and  G*  defined  in  Theorem  [8} 

( 1  - 2X2x) e“Xl  ZG 1  - ( 1  - 2Xi x) e~**z G2 -  (h -X2) xGf  =  X2-Xu 


e~XlZGi-e 


-XfflZG2  =  (X2-X1)(^G+  +  l), 


(105) 


The  problem  (|105|)  can  be  reduced  to  integral  equations  for  the  boundary  values  of  G 1  and  G2  based  on  the 
generalized  Cauchy  integral  formula. 

Theorem  9  (adjoint  integral  equations,  case  (a))  Let  M  /  0,  Rm  /  0,  and  RmR  f  M2.  In  this  case,  (|105[> 
yields  two  integral  equations  for  the  boundary  values  of  Fk(Q  =  e~XkZG, t(Q,  k  =  1,2: 


1  —  2xXk 
2x(X\  —  X2) 


(FfiQ-FfiQ)  +  WH&, x-Xu))  Fk(x)dx  =  -±  k 


=  1,2, 
(106) 


where  £  £  £  and  X]  and  X2  are  defined  in  Theorem  [2]  A  solution  to  ( 106 1  is  determined  up  to  a  real-valued 
constant  c,  i.e.  Fk(Q  =  c,  k  =  1,2,  are  a  homogeneous  solution  to  ( 106 1.  Let  Fk( Q,  k  =  1,2,  solve  ( 106 1,  then 
Gk(Q,  k  =  1,2,  are  determined  by 


where 


Gk{Q  =  (Fk{Q-c)zXkZ,  k  =  1,2,  ^<=£, 

c  =  l-Fk(Q  +  Fk(T)J*b~4  -Xk)dx,  C  G  £■ 

2  2m  Jn]l, 


(107) 


(108) 
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Proof.  The  proof  is  completely  analogous  to  the  proof  of  Theorem  [4]  □ 

In  fact,  the  integral  equations  (|5T|)  and  (|106[)  are  very  similar.  The  only  difference  is  in  the  sign  of  A*  in  the 
kernels.  The  next  proposition  shows  that  solutions  of  (|5]j)  and  (|106|)  are  closely  related. 


Proposition  11  Let  £  be  the  reflection  of  the  curve  £  over  the  r  axis,  and  let  Fk{ff),  k  =  1,2,  be  a  solution  to 
the  integral  equations  (| 5 1  [)  for  £  £  l,  then  F/ff)  =  FflC),  k  =  1 . 2,  j  £  £,  is  a  solution  to  the  integral  equations 

(fT06l>. 


Proof.  With  the  change  of  variables  x  =  X*,  X*  £  £,  and  £  =  C,*,  C,*  £  t,  in  (|106[)  and  with  the  properties 


Tf|(^,x)  =  ‘W'(^x)  and  <Hj/(£,x,  —A)  =  A),  the  integral  equations  (|106[)  reduce  to  (|5TT)  for  unknown 

functions  Fk{Zf)  =  Fk(t*),  k  =  1,2,  Cf  el.  □ 


6.2  Case  (b):  M^O  and  Rm  =  0 

We  recall  that  for  Rm  =  0,  the  magnetic  field  is  uncoupled  from  the  velocity  field  but  is  involved  in  determining 
the  velocity  field,  and  Corollary  [T]  shows  that  the  magnetic  field  is  zero  in  this  case.  For  the  adjoint  equations 
([94])— (|95|) ,  when  Rm  =  0,  w  and  q  are  uncoupled  from  g+  and  v  ,  and  their  representation  follows  directly  from 
(flQ4|). 


Corollary  4  (solution  representation  for  the  adjoint  equations,  case  (b))  Let  Ai  and  A2  be  as  those  defined 
in  Corollary  [7]  In  the  axially  symmetric  case  with  M  f  0  and  Rm  =  0,  a  solution  to  the  adjoint  equations 
is  given  by 


w-  T  iwr  =  — 


1 


Ai  —  A2 


G 1  -e~XlZG2 


q  +  ia*  =  e-x'z[  Gi  - 


R  _ 


2(Aj  —  A2) 


Gi  +e 


— X2  z 


G2  + 


R  _ 
G2 


2(Ai  —  A2) 


(109) 


v  +i  g  = 


Ai  —  A2 
v+  +  ig+  =  G+ 


A2  e-^1  z  Gi  —  Ai  q~^2ZG2]  +  Gj, 


where  G\,  G2,  and  Gf  are  defined  in  Theorem [s] 


In  this  case,  both  equations  in  (|105[)  reduce  to  the  same  boundary-value  problem 

>-ki  z, 


e  '“llG\  —  G2  =  A2  —  Ai ,  ^ei. 


(110) 


Once  Gj  and  G2  are  found  from  (|1 10[),  the  boundary  values  of  Gf  can  be  found  based  on  the  last  two  formulas 
in  ( 109 1  and  the  fact  that  v+  +  ig+  =  v~  +  ig”  on  l.  Namely,  G^  —  Gj  =  {f2e~XlZG\  —  Ai  e,~XlZG2) 


on 


£,  and  G3  are  determined  by  the  Cauchy  integral  formula  ([7])  for  r-analytic  functions. 

Theorem  10  (adjoint  integral  equations,  case  (b))  Let  M  f  0  and  Rm  =  0,  and  let  F\  (Q  =  e  ~XlZG\  (Q,  where 
Ai  is  defined  in  Corollary^  then  ( 1 10 1  reduces  to  the  integral  equation  for  the  boundary  value  of  F\ : 


1  I  feX2^WH^T,-L2)-e^z  ^WH(^T,-ll))Fl(T)dz  =  l2-ll,  (111) 


A  solution  to  (64 1  is  determined  up  to  a  real-valued  constant  c.  Let  F\  (Q  solve  (HD-  then  Gi(Q  in  ( 1 10[)  is 
given  by 

G1(Q  =  (f1(  Q-c)eH  ?£*,  (112) 
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where 


(113) 


c=lM  Q  +  ^-<P  ?€*. 
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Proof.  The  proof  is  analogous  to  the  proof  of  Theorem  [5]  □ 

Proposition  12  Let  I  be  the  reflection  of  the  curve  i  over  the  r  axis,  and  let  F\(f)  be  a  solution  to  the  integral 
equation  (|64[)  for  ?  E  I,  then  F\  (?)  =  F\  (?),  ?  E  l,  is  a  solution  to  the  integral  equation  (|1 1 1[). 

Proof.  The  proof  is  similar  to  that  of  Proposition  [IT]  □ 


6.3  Case  (c):  M^O  and  RmR  =  M2 

Theorem  11  (solution  representation  for  the  adjoint  equations,  case  (c))  Let  L  =  (R  +  Rm)/2.  In  the  axi¬ 
ally  symmetric  case  with  M  f  0  and  RmR  =  M2,  a  solution  to  the  adjoint  equations  (|94|)-(|95|)  is  given  by 


wz  +  \wr  = 


q  = 


1  ( 

2R 

R+  Rm  V 

R  +  Rm 

Rm  R 

(  2  e 

R  +  Rm 

V  R  +  Rm 

Gj  ■ 

RmR  / 

2  e-» 

-Xz  , 


-Xz , 


1  \  R 

2  /  +  R  +  R 

RTrJ^3'’ 


G2  +  RmG 


3  ’ 


(114) 


e'XlRe[Gi]+Re 


R  +  Rm  V  R  +  Rn 


1  1 

Z - r - 

2  R  +  Rn 


G2 


+  Re[G^]  +  Re[G2], 


to  = 


1 


R+Rn 


2Re~X'Tm[Gi]  +  Rm  Im[G: 


where  G\  is  an  H -analytic  function  in  D  satisfying  ([3])  with  L  and  vanishing  at  infinity;  Cl2  and  G3  are 
r-analytic  functions  in  D  ~  and  vanishing  at  infinity;  and  G,  is  an  r- analytic  function  in  D+. 


Proof.  The  proof  is  similar  to  the  proof  of  Theorem  [3] 

In  this  case,  the  boundary  conditions  (|96|)  reduce  to  the  boundary-value  problem 

*-Xz.~  ~  -  '  1  ^ 


□ 


2e-^G!-G2  +  (R  +  Rm)  ( -G+  +  1  )  =0,  ?e£, 


Z--r)G2  +  Gj-l)+G+  =  0,  ?EC 


(115) 


Theorem  12  (adjoint  integral  equations,  case  (c))  Let  M  f  0  and  Rm  R  =  M2,  and  let  L  =  (R  +  Rm)/2,  then 
the  boundary-value  problem  (|1 15  >  reduces  to  two  integral  equations  for  F\(Q  =  e  k:Gi(Q  and  F2(Q  =  G2(Q: 


_Li  (e*(«-‘)  WH(^ x-X)  -  ^(?,t))  F1(x)dx  +  Fl(Q  -  l-F2(Q  =  0, 


27ti 


tui' 


R 


MuAZ,-Z+-2(r-r'T2{Z)W^Z)dT-Rm(  R  +  Rm) 


(2  Fl(Q-F2(Q)  = 


R  +  Rn 
Rm 


(116) 


where  ?  E  £  Equations  (JTT6  >  determine  F\  and  F2  up  to  c  and  2c,  where  c  is  a  real-valued  constant.  Let  F\  and 
F2  solve  <[TT6l).  then  the  solution  to  the  boundary -value  problem  d  1 15|)  is  given  by 


G\ (?)  =  ( Fi (?)  -  c )  e^,  G2(?)  =  F2(?)  -  2c,  ?  E  l, 


(117) 
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where 


(US) 


2c=\F2(Q  +  ^-{f>  F2(z)Wr(Z,z)d%,  £e£. 
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Proof.  The  proof  is  similar  to  that  of  Theorem  [6]  □ 

Proposition  13  Let  i  be  the  reflection  of  the  curve  £  over  the  r  axis,  and  let  Fk(lf),  k=  1,2,  be  a  solution  to  the 
integral  equations  <] 7 0 [)  for  L,  E  t,  then  Fk(Q  =  Fk( Q,  k  =  1,2,  ^  E  £,  are  a  solution  to  the  integral  equations 

(TT6l>. 

Proof.  The  proof  is  similar  to  that  of  Proposition  [IT]  □ 


7  Asymptotic  Behavior  at  Conic  Endpoint 

It  is  known  that  for  a  nonconducting  viscous  incompressible  fluid  under  the  Stokes  and  Oseen  approximations, 
the  minimum-drag  shapes  subject  to  a  volume  constraint  have  conic  endpoints  with  the  angle  between  the  axis 
of  revolution  and  the  tangent  at  the  endpoint  to  be  27t/3;  see  j20l  [171.  This  section  analyzes  the  asymptotic 
behavior  of  a  solution  to  the  MHD  problem  <fT0|)— <fT3|)  near  the  vicinity  of  a  conic  endpoint. 

Let  one  of  the  endpoints  of  £  lie  on  the  z  axis  and  have  coordinates  r  =  0,  z  =  c,  and  let  (p,  0)  be  the  local 
polar  coordinates  with  the  pole  at  r  =  0,  z  =  c  and  with  the  angle  0  counted  from  the  z  axis.  In  this  case,  (p,  0) 
are  related  to  (r,z)  by  r  =  psin0  and  z  =  c  +  pcos0,  0  E  [0, 7c] .  Let  0o  be  the  angle  between  the  z  axis  and  the 
tangent  to  £  at  the  endpoint.  It  is  assumed  that  £  is  in  the  cone  0  >  0q. 


7.1  Cases  (a)  and  (b):  M^O  and  Rm  R  M2 

Proposition  14  (asymptotic  behavior  of  G\  and  Gf)  Let  M  f  0  and  Rm  R  f  M2,  and  let  G\  and  Gi  be  H- 

analytic  functions  defined  in  Theorem  [2]  The  asymptotic  behavior  of  G\  and  G2  on  t  in  the  vicinity  of  the 
endpoint  is  given  by 

(  bk  +  o( p°),  0O  <  2tc/3, 

Gk=  l  ak  (lnp  +  21n(cos(0o/2))  —  i  tan(0o/2) )  +  bk  +  o  (p°)  ,  0o  =  27t/3,  (119) 

[  flr-p“~1Ma_i(cos0o)  +  ^  +  o(po) ,  0O  >  2ji/3, 


as  p  — >0,  where  ak,bk  El  are  some  constants,  Ma(cos0)  =  (a+  I )  Pa(cos0)  +iP«  '(cos0),  and  (X  is  the  single 
zero  of  the  equation 


[  I  +  cos2  0o )  Pa  (cos  0o )  Pa  ^  (cos  0o )  + sin  0o  cos  0o 


Pa*  (cos  0O) 


i  2 


+  a(a+l)[Pa(cos0o)f  =0  (120) 


in  the  interval  (0,  \  )  for  given 


Proof.  The  proof  is  conducted  similarly  to  the  proof  of  Proposition  1  in  1(321 .  First,  let  Rm  f  0. 

In  the  local  polar  coordinates  (p , 0),  the  //-analytic  functions  G\  and  G2  and  r-analytic  functions  G^  can 
be  represented  in  the  form 

Gk  =  arpa_1Ma„i(cos0)  +bk  +  pa  (c*Ma(cos0)  —  /,^fl^Pa(cos0))  +  o(pa),  0  <  a  <  1, 

G±  =  p^'Ma-^cos©)  +  bf  +  c^paMa(cos0)  +  o(p“),  0  <  a  <  1, 

7  Equation  <]  1 20[)  reduces  to  identity  for  a  =  0  and  a  =  1 . 
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as  p  — >  0,  where  k  =  1,2,  and  a\,  b\,  c i,  a2,  bi,  C2,  af,  bf,  cf  are  real  constants.  Substituting  d  1 2 1 1  into  the 
boundary-value  problem  (|38|)  with  eXkZ  =  e^*c(l  +  /+p  cos0  +  0(p2))  and  equating  corresponding  coefficients 
at  pa_1,  p°,  and  p“,  we  obtain  six  equations: 

(1  -2X2x)flie^lC-(l  -2^ix)a2e:>l2C-R'>llR_'>^|2  a~  =  0, 

(1  —  2X2^)b\&x'c  —(1  —  2X\ x)62e^2'  -  R  R — ^j2  ^3  =  ^2  —^i, 

e^1  c  a i  —  eXlC  a2  =  (fh2  —  Xi  )af , 

eh c bx  - eXl c  £2  =  (X2  - ^ ) (bf  + 1 ) , 

(1  -2X2x)Cle^c-(l  -2Xix)c2eX2C-  RXR~_^2  c3~)  M«(0 

+  ^(1  —  2X2x)Xi  fli  e^lC  —  (1  —  2Xix)X2a2e^2^  (/Mtt_i(i)  —  P«(0)  =  0; 
(e^cci-eX2Cc2-(X2-X1)4)Ma(0  +  (e^cX1fl1-eX2CX2fl2)  (tMa_!(t)  -P„(/))  =  0, 

where  t  =  cos0o.  In  fact,  (|122[)  and  ( 123[)  are  complex-valued  equations  and  equivalent  to  a  system  of  four 
equations  which  can  have  nonzero  a i  and  a2  only  if  system’s  determinant  is  zero,  i.e., 


(122) 


(123) 


Re[M a(t)]  Im  [t  Ma_!  (t)  -  Pa(t)]  -  Im[Ma(f)]  Re  [f  M^f)  -  Pa(f)]  =  0. 


(124) 


With  the  relationships  Pa-i(f)  =  fPa(f)  —  ON1  \/l  —t2p{)'(t)  and  P„2j(f)  =  /P«^(f)  +  a y/\  —  t2Va(t),  the 
condition  (|124[)  reduces  to  (|120[),  which  has  a  zero  a  in  (0, 1)  only  for  0o  >  2ji/3.  The  zero  is  single  on  (0, 1), 
and  a  — >  1—  as  0o  — >  271/3+. 

Now  when  a  — >  1,  pa_1Ma_i(cos0)  becomes  equivalent  to  p°  in  ( 121 1.  In  this  case,  G\,  G2,  and  Gf  are 
represented  by 

Gk  =  (A*(0)  +-B*(0)  p)  Inp  +  Q(0)  +Dk(d)  p  +  o(p),  *=1,2, 

G3  =  {Af(d)  +  Bf(d)p)lnp+Cf(d)+Df(d)p  +  o(p), 
with  complex-valued  functions  A\,B\,  C\,  A2,  fi2,  C2,  A*,  Bf,  Cf  given  by 


(125) 


Ak(Q)=akl  Bk(d)  =  bkMl(cosQ)-aklkcose,  Ck(d)  =  ck  +  ak\ j//(0),  *  =  1,2, 

Dk(Q)  =  dk Mi  (cos 0)  —  ^\|//(0)  +Xk  (ak  (l  —  2 cos 0  In  (cos  §))  —  ck  cos0)  ,  k  =  1,2, 

A±(0)  =  af,  Bf(8)  =  ^Mi(cos0),  C^(0)  =  cf  +a^\j/f  (0),  D±(0)  =  rf±Mi(cos0) ~bf\\rf  (0) 

where 

\|//(0)  =  21n(cos(0/2))  —  i  tan(0/2),  \|//  (0)  =  —  \|//(0)Mi(cos0)  +  1  —  icos0  tan(0/2) 

\|/j"(0)  =21n(sin(0/2))  +  icot(0/2),  \|//(0)  =  —  \|//(0)Mi(cos0)  —  1 +  icos0  cot(0/2), 

and  a\,  b\,  ci,  a2,  Z?2,  c2,  af,  bf,  cf  are  real  constants.  Observe  that  Ai,  B\,  C\,  A2,  R2,  C2,  A/,  R/,  C f  are 
finite  for  0  E  [0, 7c),  whereas  Af,  Bf,  Cf  are  finite  for  0  E  (0,Jt]. 

For  brevity,  let  ak  =  (ak,bk,ck),  k  =  1,2,  and  a/  =  (af  ,bf  ,cf).  Substituting  (|125|)  into  (|38j)  with  eXkZ  = 
e+c(i  +  X,tpcos0  +  0(p2))  and  equating  coiTesponding  coefficients  at  lnp,  plnp,  p°,  and  p,  we  obtain  after 
some  transformations:  af  =  0,  bf  =  0,  and 

(1  -  2X2x)  e^lC  ai  -  ( 1  -  2A,i x)  &XlC  a2  -  R^1R_^2  a/  =  (0, 0,  -  X2) , 


e?L1C  3]  —  eA2C  32  —  (A,2  —  =  (0,0,  X2  —  Xi), 


=+2  c 
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( 1  -  2\2k)  di  eX'c  -( 1  -  2Xi  x)  d2  el2C 


(126) 


-R^j5dOMl(COS6o) 

+  ^(1  —  2X2x)Xi  a\  e^lC— (1  —  2Xix))i2fl,2e?l2C^  (1  —  icosGo  tan(0o/2))  =  0, 

(di  eXlC  -d2  el2C  ~(M  -  h )  d£\  Mi  (cos  0O)  +  (li  a\  &X{C  -X2  a2  e^2C)  ( 1  -  i  cos  0O  tan(0o/2) )  =  0.  (127) 

As  in  the  previous  case,  (|  1 26}-([T27}  is  a  system  of  four  scalar  equations  which  can  have  nonzero  a\  and  a2 
only  if  system’s  determinant  is  zero: 

Re[Mi(cos0o)]Im[l  —  icos0o  tan(0o/2)]  —  Im[Mi(cos0o)] Re[l  —  icos0o  tan(0o/2)]  =  0,  (128) 

which  reduces  to  2cos20o  —  cos0o  —  1  =  0.  This  condition  implies  0o  =  0  and  0o  =  2tt/3,  which  proves  (|1 19[) 
for  0o  =  2ji/3.  For  Rm  =  0,  the  proof  is  completely  analogous.  □ 


7.2  Case  (c):  M  ±  0  and  Rm  R  =  M2 

Proposition  15  (asymptotic  behavior  of  Gi  and  G2)  Let  M  /  0  and  Rm  R  =  M2,  and  let  G\  and  G2  be  those 
defined  in  Theorem  [i]  The  asymptotic  behavior  of  G\  and  G2  on  t  in  the  vicinity  of  the  endpoint  is  given  by 
(fIT9l)-(fl20l). 


Proof.  The  proof  is  similar  to  that  of  Proposition  [14]  Though  there  are  some  differences.  In  the  vicinity  of 
the  conic  endpoint,  let  G\  and  G't  be  represented  by  (|121|)  with  k\  =  ( R  +  Rm)/2.  However,  G2  here  is  an 
r-analytic  function  and  is  represented  by 

G2  =  a2pa_1Ma_i(cos0)  +  &2  +  c2paMa(cos0)  +  o(pa),  0  <  a  <  1. 

Substituting  these  representations  into  (|67j)  and  equating  corresponding  coefficients  at  pa,  p°,  and  p“,  we  have 

2ni  e^l<;  —  a2  +  (R  +  Rm)  a^  =  0,  Rm  (ca2  +  af)  —  R«^  =  0, 

2b\  tx'c  —b2  +  (R  +  Rm)  (bj  +  1)  =  0,  Rm  {cb2  +  bf  +  l)  -  Rb+  =  0, 


2 ci  eX|C  -c2  +  (R  +  Rm)  c+)  Ma(?)  +  2\iai  eXlC  (? Ma_i  (?)  -  Pa(0)  =  0, 
(Rm  {cc2  +  cf)  -  Rc^)Ma(f)  + ^Rm<32M1(f)Mct_i(t)  =  0, 


(129) 


where  t  =  cos0o-  The  complex- valued  equations  (|129[)  have  nonzero  a\  and  a2  if  the  determinant  of  each 
equation  is  zero.  The  determinant  of  the  first  equation  in  d  129[)  is  given  by  d  1 24|),  whereas  the  determinant  of 
the  second  takes  the  form 


Re[Ma(t)]  Im  [Mi (?)  Ma_i  (?)]  -  Im[Ma(?)]  Re  [Mi  (?)  Ma_i  (?)]  =  0. 


Remarkably,  as  (124},  this  condition  also  reduces  to  (|120|),  which  proves  (|119)  for  0  <  a  <  1. 

When  a  — r  1,  the  functions  Gi,  G2,  and  Gf  are  represented  by  (|125[)  with  Xi  =  (R+  Rm)/2  and  with 
complex- valued  functions  Ai,  B\,  C\,  Bf,  as  defined  in  the  proof  of  Proposition  [l~4|  and 

A2(Q)  =  a2,  B2(B)  =  b2M{(cosQ),  C2(0)  =  c2-Fa2Vr(®)’  D2(B)  =  d2Mi(cos0)  —  b2yf  (0). 


Substituting  these  representations  into  ([67}  and  equating  corresponding  coefficients  at  lnp,  plnp,  p°,  and  p,  we 
obtain  after  some  transformations:  a ^  =  0,  bf  =  0,  2e^lCfli  —  a2  =  0,  2eXlCb\  —  b2  =  0,  2e^lCci  —  c2  +  (R  + 
Rm)(c3~  +  1)  =  0,  ca2  +af  =  0,  a2/2  +  cb2  +  bf  =  0,  Rm  (cc2  +  cf  +  l)  —  Rc^  =  0,  and 


^2e^lCdi  —d2^j  Mi  (cos  0o)  +2^i  a\ &XlC  (1  —  i  cos0o  tan(0o/2))  =  0, 
(Rm(c2/2-| -cd2  +  df)  —  Rd^)  Mi  (cos  0o )  +  (1  —  icos0Q  tan(0o/2))  =  0. 


(130) 
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The  complex-valued  equadons  ( 130  >  have  nonzero  ai  and  a 2  if  ( 128 1  holds,  and  the  rest  of  the  proof  is  similar 
to  that  of  Proposition  [14]  □ 


7.3  Adjoint  Equations 

The  asymptotic  behavior  of  solutions  to  the  adjoint  equations  ([94])— ([96])  near  the  conic  endpoint  on  the  boundary 
is  analyzed  similarly. 

Proposition  16  (asymptotic  behavior  of  G\  and  Gj  for  the  adjoint  equations)  Let  M  f  0,  and  let  G\  and  Gn 

be  defined  as  in  either  Theorem^or  Theorem [77]  In  both  cases,  the  asymptotic  behavior  of  G \  and  Gi  on  I  in 
the  vicinity  of  the  endpoint  is  given  by  (fTT9l)-(fT20|). 

Proof.  The  proof  is  completely  analogous  to  the  proofs  of  Propositions  [l4|  and  [15]  □ 


8  Analysis  of  Minimum-drag  Shape 


The  results  obtained  in  Sections[5]-[7]have  the  following  implications. 

With  the  representations  (|104[),  (|109[),  and  (|1 14[),  we  can  prove  that  in  PropositionfTO]  Ir  — >  0  as  R  — >  °o.  The 
proof  is  similar  to  showing  that  Ir  — r  0  as  R  — >  °o  in  the  proof  of  Proposition^  and  consequently,  the  assumption 
of  Proposition[TOlholds  true.  With  the  representations  for  solutions  of  the  adjoint  equations  obtained  in  Section 
[6]  the  optimality  condition  (|103[)  is  specialized  for  each  case  (a),  (b),  and  (c). 


Corollary  5  (optimality  condition) 

(a)  In  case  (a)  (i.e.  M  f  0,  Rm  R  f  M2,  Rm  f  0),  let  G\  and  G7  be  those  defined  in  Theorem [2]  and  let  G\  and 
Gj  be  those  introduced  in  Theorem^  The  representations  (|18[)  and  ( 1Q4[)  with  the  boundary  conditions 
(38 )  and  ( 105 )  and  with  the  identity  (R  —  2A.i)(R-  Tkf)  +  M2  =  0  imply  that  on  i,  the  optimality  condition 
( 103|)  simplifies  to 


coco*  +  Rm/rr  g  = 


Xl  —  h 


((R  —  I'Ll)  Im[Gi]  Im[G*]  -  (R  -  2X1)Im[G2]  InifG,])  =  const.  (131) 


(b)  In  case  (b)  (i.e.  M  f  0,  Rm  =  0),  let  G\  and  G2  be  those  defined  in  Corollary  [7J  and  let  G\  and  G\  be 
those  introduced  in  Corollary  [7]  The  representations  (|29[)  and  ( 109[)  with  the  boundary  conditions  (60 ) 
and  ( 1 10|)  imply  that  on  i,  the  optimality  condition  (|103[)  reduces  to 


coco*  +  Rm/ir  g  =  4Inr[Gi]  InrfGj]  =  const. 


(132) 


and  Im[Gj]  =  ImfGj]  on  l,  and  under  these  conditions,  ( 1 32 1 


Observe  that  in  this  case,  Inr[Gi]  =  Im[G2] 
follows  from  (|  1 3 1 1). 

(c)  In  case  (c)  (i.e.  M  f  0,  Rm  R  =  M2),  let  G\  and  G2  be  those  defined  in  Theorem^  and  let  Gj  and  Gj  be 
those  introduced  in  Theorem  11  The  representations  ©  and  (|1 14|)  with  the  boundary  conditions  (|67[) 
and  ( 1 15[>  imply  that  on  i,  the  optimality  condition  ( 103[)  simplifies  to 


coco*  +  R  mhrg  = 


1 

R  +  Rn 


(4R  Im[Gi]  Im[Gj]  +  Rm  Inr[G2]  Im[G|])  =  const. 


(133) 


Corollary  6  If  Rm  =  0  and  R  =  0,  the  adjoint  equations  (|94| )  for  w  and  q  coincide  with  (|13[)  for  u  and  p.  In  this 
case,  the  boundary  conditions  (|10[)  and  \96\  imply  w  =  — u  and  CO*  =  —CO,  and  thus,  the  optimality  condition 
( 103|)  simplifies  to  0)  =  const  on  S,  or  eMz/2  Im[Gi]  =  const  on  t,  where  G\  is  defined  in  Corollary  [2] 
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Corollary  7  Propositions  \  1  7f[7j|  imply  that  if  £  is  fore-and-aft  symmetric,  i.e.  symmetric  with  respect  to  the 
r  axis,  and  l  admits  a  parametrization  C,  =  l '(t),  t  E  [—1, 1],  such  that  £(f)  =  £(— t),  then  if  Ff,(t),  t  E  [—1, 1], 


k  =  1,2,  are  a  solution  to  (|51[),  or  ([64]),  or  m  then  Fk(—t),  t  E  [—1, 1],  k  =  1,2,  are  a  solution  to  d  106[)  or 
(|  1 1 1 1)  or  (11 6|),  respectively. 

Now  if  an  optimal  shape  is  to  be  found  iteratively  and  the  iteration  process  starts  from  a  fore-and-aft 
symmetric  shape,  e.g.  sphere  or  spheroid,  then  at  each  iteration,  the  optimality  conditions  (|1 3 1 1— ( 133 1  are 
symmetric  with  respect  to  the  r  axis,  which  means  that  the  optimal  shape  will  be  fore-and-aft  symmetric.  In  this 


case,  ( 131 1— ( 133 1  take  the  form  l)(t)  =  const  for  t  E  [— 1, 1],  where 


v(t)  =  < 


A.1  —  )i2 
-4Vi(/)Vi(-/) 

1 


R  +  Rn 


((R  —  2A,2)Vi(f)  Vi(— t)  —  (R  —  2'ki)V2(t)  V2(— t))  in  case  (a), 

in  case  ( b ), 

(4R V\(t)  V\ (— f)  +  Rm  V2 (t)  V2 (— f))  in  case  (c), 


(134) 


and  14(f)  =  \m[Fk(t)\,  k=  1,2. 


Propositions  [T4|  and  15  do  not  specify  whether  a\  and  a 2  for  Gi  and  Gi  in  (]  1 19|)  arc  related.  Since  the 
kernels  in  the  integral  equations  (|5T|)  and  ([70])  have  logarithmic  singularities,  the  integrals  in  ([51])  and  (]70|)  with 
G'i  and  G?  behaving  as  (JTT9J)  in  the  vicinity  of  a  conic  endpoint  are  integrable  and  finite.  This  implies  that 
F]  (Q  —F2(Q  in  ([51])  and  2F\  (Q  —  G2&)  in  (f70[)  are  finite,  and  consequently,  in  (|119),  ai  =  02  for  cases  (a)  and 
(b),  and  2a  1  =  a 2  for  case  (c).  The  same  is  true  for  Gj  and  Gf 

Corollary  8  (conic  endpoints)  In  the  setting  of  Corollary^  suppose  t  =  ±1  correspond  to  conic  endpoints. 
Then  in  the  vicinity  oft  =  ±1,  v>(f)  behaves  as  — 414(f)  Vi(— f)  in  all  three  cases  (a)-(c)  with  V\  being  defined 
respectively.  Then  ( 1 19|)  implies  that  —  414(f)  V\  (— f)  =  const  f  0  in  the  vicinity  o/±l  only  iff}  0  =  2ft/3. 


Now  we  proceed  with  finding  the  minimum-drag  shape  in  semi-analytical  form.  Let  C,  be  parametrized  by 

?(0  =  K0+iz(0  =y(cos(ro/2)£"lo«; T2j{t)  +  ism(m/2)Yjnj1=0bjT2j(t)^  ,  t  E  [-1,1],  (135) 

where  Tj(t)  is  the  Chebyshev  polynomial  of  the  first  kind;  ao,...,ani,  bo,...,bn2  are  real  coefficients;  and 
y  =  y(ao,  ■  •  •  ,ani,bo, . . .  ,b„2)  is  the  multiplier  introduced  to  satisfy  the  volume  constraint  k  f_1  r2z!dt  =  4ji/3 
identically,  i.e.  if  Cft)  =  y (?(f)  +iz(f)),  then 


Y  = 


'r'zdt 


-1/3 


(136) 


The  choice  of  the  parametrization  (|  1 35  )-( 136 )  is  explained  in 

Let  p(f)  =  'u(f)  —  j  J_*j  v(t)dt,  where  the  second  term  in  the  right-hand  side  is  the  average  of  v(t)  on 
f  E  [—1,1],  and  let  ^  =  4  +  i Zk  and  0/c  be  the  shape  ^  and  function  o,  respectively,  at  step  k.  Then  i can  be 
updated  as  in  Pironneau’s  iterative  procedure  |i20ll : 


^k+l  =  t>k  +  £kV>k 


(137) 


This  procedure  reduces  to  finding  optimal  ao,-.-,ani,  bo,.. .  ,b„2  in  (|135|)  as  described  by  Algorithm  1  in 
031.  in  which  u*  needs  to  be  replaced  by  u,  and  initial  «o  and  bo  should  correspond  to  the  axes  of  the  minimum- 
drag  spheroid  found  in  Section  |4~2]  (for  given  R,  Rm,  and  S)  while  initial  ai, . . .  ,aHl,  b\.. . . .  b„2  should  be  set 
zero. 
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Figure  [4ji  shows  three  minimum-drag  shapes,  having  the  volume  of  the  unit  sphere,  for  R  =  Rm  =  3  and 
S  =  0,  1,  2.  In  the  case  of  S  =  0,  which  corresponds  to  the  Oseen  equations,  the  minimum-drag  shape  was 
obtained  by  the  approach  presented  in  lf32l.  For  the  three  shapes,  Tables  [T]  and  [2]  present  corresponding  drag 
coefficient  Cd,  the  constant  in  the  optimality  condition  (|103[),  error  ||  -  ||  in  satisfying  (|103|),  and  parameters  in 
the  representation  (|135[)-(fT36l).  Among  those  three  shapes,  the  one  for  S  =  1  is  shortest.  Let  C*D  be  the  drag 
coefficient  for  the  unit  sphere.  Figure  [4|r  interpolates  the  ratio  Cd/C*d  for  R  =  Rm  =  3  based  on  three  obtained 
values  for  S  =  0, 1,2.  Figures^  and[4p  show  that  the  drag  ratios  CujC*D  for  the  minimum-drag  spheroids  and 
minimum-drag  shapes  for  same  R,  Rm,  and  S  are  sufficiently  close. 


-2-10  1  2 


(a)  cross  section  of  minimum-drag  shapes 


Figure  4:  Minimum-drag  shapes  and  the  interpolated  drag  ratio  Cd/Cq  for  R  =  Rm  =  3  and  S  =  0, 1,2. 


Table  1 :  Drag  coefficient  Cd,  the  constant  in  the  optimality  condition  (|103|),  and  error 
for  the  minimum-drag  shapes  for  R  =  Rm  =  3  and  S  =  0,  1,  and  2. 


in  satisfying  (|103|> 


S  =  0 

S=  1 

S  =  2 

CD 

1.64828 

1.36736 

1.81103 

constant 

-3.63472  x  KT1 

-2.35380 

-3.66957 

error 

1.58348  x  1(T3 

3.94289  x 10^3 

4.91428  x  10~3 

9  Conclusions 

This  work  has  developed  an  approach  of  generalized  analytic  functions  to  the  MHD  problem  of  an  electrically 
conducting  viscous  incompressible  flow  past  a  solid  body  of  revolution  under  the  assumptions  that  the  applied 
magnetic  field  and  body’s  axis  of  revolution  are  aligned  with  the  flow  at  infinity  and  that  the  body  and  fluid  have 
same  magnetic  permeability.  In  the  three  complementary  cases:  (a)  Rm  /  0,  RmR  /  M2  (S  /  1);  (b)  Rm  =  0; 
and  (c)  RmR  =  M2  (S  =  1)  all  assuming  M  /  0,  the  velocity  field,  pressure  and  magnetic  fields  in  and  out  the 
immersed  body  are  represented  by  four  generalized  analytic  functions.  In  particular,  in  case  (a),  the  fields  are 
represented  by  two  //-analytic  functions  and  two  r-analytic  functions;  in  case  (b),  the  magnetic  field  is  constant 
everywhere,  and  the  velocity  and  pressure  are  represented  by  two  //-analytic  functions;  and  in  case  (c),  the 
fields  are  represented  by  one  //-analytic  function  and  three  r-analytic  functions.  r-Analytic  functions,  whose 
real  and  imaginary  parts  are  harmonic,  are  related  to  the  Stokes  equations,  whereas  //-analytic  functions  have 
the  real  and  imaginary  parts  satisfying  the  modified  Helmholtz  equation  and  are  related  to  the  Oseen  equations. 
Namely  the  difference  in  the  number  of  r-analytic  and  //-analytic  functions  in  the  solutions  in  cases  (a)  and  (c) 
contributes  to  the  peculiarity  of  the  case  of  S  =  1 . 

In  all  three  cases  (a),  (b),  and  (c),  the  MHD  problem  has  been  reduced  to  integral  equations  for  unknown 
boundary  values  of  the  involved  generalized  analytic  functions  based  on  the  generalized  Cauchy  integral  for- 
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Table  2:  Parameters  in  the  representation  ( 135 )-( 136 1  for  the  minimum-drag  shapes  for  R  =  Rm  =  3  and  S  =  0, 
1,  and  2. 


o 

II 

S  =  1 

S  =  2 

Y 

1.01196396246 

1.01366399904 

1.01131424639 

«o 

6.29221950407  x lO^1 

6.50248611262  x  10_1 

5.93267526768 x 10“ 1 

a\ 

-9.43328552499  x  10~2 

-9.23964575126 x 10-2 

-1.17578599834  x  10'1 

Cl2 

-1.37355103545  x  10~2 

-6.7092264183  x  10~3 

-2.61378349707  x  10~2 

«3 

1.98202204631 x 10~3 

2.34809221062  x  10~3 

-1.57412245298  x  10~3 

CI4 

- 

4.22269308778  x  10~4 

-3.89459211608 x  10~5 

- 

- 1.08999300292 x  10'4 

-2.62129253651 x  10~5 

bo 

2.00004216653 

1.87568864891 

2.13349213677 

b\ 

8.69011301008 x 10~2 

1.00707086887  x 10_1 

8.57434199435  x 10~2 

b2 

1.39653034146 x 10-4 

2.70040905434  x  10~4 

-1.52613634178  x  10'3 

b-i 

2.83047975928  x  10~3 

1.85753309456 x 10~3 

1.31075784337  x 10~3 

b4 

- 

-1.0430097082  x  10~4 

-1.10494890615  x  10~5 

b  5 

- 

-7.96297414857  x 10~5 

-8.91048112023  x  10~6 

mula.  Solutions  to  the  integral  equations  have  been  represented  by  a  finite  functional  series  with  series  co¬ 
efficients  determined  by  quadratic  error  minimization  and  have  been  shown  to  coincide  with  the  series-form 
solutions  for  sphere  with  high  accuracy. 

The  necessary  optimality  condition  ( 103[)  for  minimum-drag  shapes  subject  to  the  volume  constraint  has 
been  obtained  by  Mironov’s  shape  variation  approach.  In  the  absence  of  the  magnetic  field,  it  reduces  to 
those  for  the  Stokes  and  Oseen  equations,  and  coincides  with  the  one  for  the  Navier-Stokes  equations.  It  has 
been  shown  analytically  that  regardless  of  magnitudes  of  R,  Rm,  and  S,  the  minimum-drag  shapes  have  conic 
endpoints  with  the  angle  of  27t/3  counted  from  the  axis  of  revolution.  Remarkably,  the  minimum-drag  shapes 
for  the  Stokes,  Oseen,  and  Navier-Stokes  equations  have  the  same  property.  The  numerical  analysis  has  shown 
that  the  minimum-drag  shapes  are  fore-and-aft  symmetric  and  that  for  R  =  Rm  =  3,  the  minimum-drag  shape 
for  S  =  1  has  the  largest  drag  ratio  and  smallest  aspect  ratio.  Also,  for  same  R,  Rm,  and  S,  the  drag  coefficients 
for  the  minimum-drag  spheroids  and  minimum-drag  shapes  are  sufficiently  close. 
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